Nonzero orbital angular momentum superfluidity in ultracold Fermi gases 



o 
o 

(N 



O 

o 
u 



X3 

o 



> 

O 

o 



I 

o 



M. Iskin and C. A. R. Sa de Melo 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, 

(Dated: February 6, 2008) 



USA 



We analyze the evolution of superfluidity for nonzero orbital angular momentum channels from 
the Bardeen-Cooper-Schrieffer (BCS) to the Bose-Einstein condensation (BEC) limit in three di- 
mensions. First, we analyze the low energy scattering properties of flnite range interactions for 
all possible angular momentum channels. Second, we discuss ground state (T — 0) superfluid 
properties including the order parameter, chemical potential, quasiparticle excitation spectrum, 
momentum distribution, atomic compressibility, ground state energy and low energy collective exci- 
tations. We show that a quantum phase transition occurs for nonzero angular momentum pairing, 
unlike the s-wave case where the BCS to BEC evolution is just a crossover. Third, we present a 
gaussian fluctuation theory near the critical temperature (T = Tc), and we analyze the number of 
bound, scattering and unbound fermions as well as the chemical potential. Finally, we derive the 
time- dependent Ginzburg-Landau functional near Tc, and compare the Ginzburg-Landau coherence 
length with the zero temperature average Cooper pair size. 



PACS numbers: 03.75.Ss, 03.75.Hh, 74.25. Bt, 74.25. Dv 



I. INTRODUCTION 

Experimental advances involving atomic Fermi gases 
enabled the control of interactions between atoms 
in different hyperfine states by using Feshbach reso- 
nancesiiSiiiiSiSii. These resonances can be tuned via 
an external magnetic field and allow the study of dilute 
many body systems with fixed density, but varying inter- 
action strength characterized by the scattering parame- 
ter tti. This technique allows for the study of new phases 
of strongly interacting fermions. For instance, the recent 
experiments from the MIT groupS. marked the first obser- 
vation of vortices in atomic Fermi gases corresponding to 
a strong signature of superfluidity in the s-wave {£ = 0) 
channel. These studies combinediiSiSiiiSiSiLS correspond 
to the experimental realization of the theoretically pro- 
posed Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein 
condensation (BEC) crossoveriSiifliiiiiSii^ in three dimen- 
sional (3D) s-wave superfluids. Recent extensions of 
these ideas include trapped fermionsi^ti^ and fermion- 
boson modcls^^iiii^. 

Arguably one of the next frontiers of exploration in 
ultracold Fermi systems is the search for superfluidity 
in higher angular momentum states {i ^ 0). Substantial 
experimental progress has been made recentljii^iSSiSiiSSiS^ 
in connection to p-wave {i = 1) cold Fermi gases, making 
them ideal candidates for the observation of novel triplet 
superfluid phases. These phases may be present not only 
in atomic, but also in nuclear (pairing in nuclei), astro- 
physics (neutron stars), and condensed matter (organic 
superconductors) systems. 

The tuning of p-wave interactions in ultracold Fermi 
gases was initially explored via p-wave Feshbach reso- 
nances in trap geometries for in Kefi^^ and ^Li 
in Refi2L2^. Finding and sweeping through these res- 
onances is difficult since they are much narrower than 
the s-wave case, because atoms interacting via higher 
angular momentum channels have to tunnel through a 



centrifugal barrier to couple to the bound state2£. Fur- 
thermore, while losses due to two body dipolariSLSi or 
three-bodjiiSiSS processes challenged earlier p-wave ex- 
periments, these losses were still present but were less 
dramatic in the very recent optical lattice experiment in- 
volving '"^K and p-wave Feshbach resonances^^. 

Furthermore, due to the magnetic dipole-dipole in- 
teraction between valence electrons of alkali atoms, the 
nonzero angular momentum Feshbach resonances corre- 
sponding to projections of angular momentum £ [mi = 
zt£,zL{£ — 1 ),..., 0] are nondegenerate (separated from 
each other) with total number of £ -I- 1 resonances'^^. 
Therefore, in principle, these resonances can be tuned 
and studied independently if the separation between 
them is larger than the experimental resolution. Since 
the ground state is highly dependent on the separation 
and detuning of these resonances, it is possible that p- 
wave superfluid phases can be studied from the BCS to 
the BEC regime. For sufficiently large splittings, it has 
been proposedi2S*2& that pairing occurs only in = 
and does not occur in the mi = ±1 states. However, for 
small splittings, pairing occurs via a linear combination 
of the mg — and mi = ±1 states. Thus, the mi = 
or mi — ±1 resonances may be tuned and studied inde- 
pendently if the splitting is large enough in comparison 
to the experimental resolution. 

The BCS to BEC evolution of d-wave {£ = 2) super- 
fluidity was discussed previously in the literature using 
cont inuum'^ ^ ' ^ ^ ' ™ and lattice^Ml descriptions in connec- 
tion to high-Tc superconductivity. More recently, p-wave 
superfluidity was analyzed at T = for two hyperflne 
state (THS) systems in SD^^^, and for single hyperflne 
state (SHS) systems in two dimensions (2D)^24^^i2^, us- 
ing fermion-only models. Furthermore, fermion-boson 
models were proposed to describe p-wavc superfluidity 
at zerc--'^^ and finite temperature-^ in 3D. 

In this manuscript, we present a generalization of the 
zero and finite temperature analysis of both THS pseudo- 
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spin singlet and SHS pseudo-spin triplet™ superfluidity 
in 3D within a fermion-only description. Our main re- 
sults are as follows. 

Through an analysis of the low energy scattering am- 
plitude within a T-matrix approach, we find that bound 
states occur only when the scattering parameter ag > 
for any £. The energy of the bound states Eh^e involves 
only the scattering length uq for £ = 0. However, an- 
other parameter re related to the interaction range 1/fco 
is necessary to characterize Et,/ for £ 7^ 0. Therefore, all 
superfluid properties for £ ^ depend strongly on ko and 
ae, while for £ — they depend strongly only on gq but 
weakly on k^. 

At zero temperature (T — 0), we study the possibil- 
ity of a topological quantum phase transition in £ ^ 
atomic Fermi gases during the evolution from BCS to 
BEG regimo^5,27,28,33,34,35,37,38^ gj^o^ ^^^^ ^^leie is 

a fundamental difference between the £ = and £ 7^ 
cases. In the s-wave {£ = 0) case, there is no phase 
transition as the magnetic field is tuned through the Fes- 
hbach resonance from the BCS to the BEC limit. That 
is, the zero temperature thermodynamic properties are 
analytic functions of the scattering length ao when the 
Feshbach resonance is crossed. In this case, the super- 
fluid ground state does not change in any fundamental 
way as ao is varied. This has been noted in the condensed 
matter literature long agoiSiiiiiSii^ and it is referred to 
as the BCS-BEC crossover problem. However, for £ ^ 0, 
we show that there is a phase transition as the magnetic 
field is swept through the £ Feshbach resonance. The 
phase transition does not occur when two body bound 
states are first formed, but occurs when the many body 
chemical potential crosses a critical value. 

To show that such a zero temperature (quantum) phase 
transition occurs in £ ^ 0, we calculate the order param- 
eter, chemical potential, quasiparticle excitation spec- 
trum, momentum distribution, atomic compressibility, 
low energy collective excitations and average Cooper pair 
size as a function of ae, and show that they are non- 
analytic at T = when the chemical potential iie crosses 
a critical value. The symmetry of the order parameter re- 
mains unchanged through the transition, as the ground 
state wavefunction experiences a major rearrangement 
of its analytic structure. In addition, the elementary ex- 
citations of the superfluid also change from gapless in 
the BCS side to fully gapped in the BEC side leading 
to qualitatively different thermodynamic properties in 
both sides. Thus, we conclude that there is a potentially 
observable BCS-BEC phase transition in £ atomic 
Fermi gases in contrast to the BCS-BEC crossover al- 
ready found in s-wave {£ = 0) gases. 

At finite temperatures, we develop a gaussian fluctu- 
ation theory near the critical temperature (T « Tc.e) to 
analyze the number of unbound, scattering and bound 
fermions as well as the chemical potential. We show that 
while the saddle point number equation is sufficient in 
weak coupling where all fermions are unbound, the fluc- 
tuation contributions have to be taken into account in or- 



der to recover the BEC physics in strong couplings where 
all fermions are bound. 

We also derive the time-dependent Ginzburg-Landau 
(TDGL) functional near Tcj and extract the Ginzburg- 
Landau (GL) coherence length and time. We recover the 
usual TDGL equation for BCS superfluids in weak cou- 
pling, whereas in strong coupling we recover the Gross- 
Pitaevskii (CP) equation for a weakly interacting dilute 
Bose gas. The TDGL equation exhibits anisotropic co- 
herence lengths for £ 7^ which become isotropic only in 
the BEC limit, in sharp contrast to the £ = case, where 
the coherence length is isotropic for all couplings. Fur- 
thermore, for any £, the GL time is a complex number 
with a larger imaginary component for ^e > reflecting 
the decay of Cooper pairs into the two-particle continuum 
with short lifetimes. However, the imaginary component 
vanishes for < and Cooper pairs become stable with 
long lifetimes above Tc/. 

The rest of the paper is organized as follows. In Sec.lTll 
we analyze the interaction potential in both real and mo- 
mentum space for nonzero orbital momentum channels. 
We introduce the imaginary-time functional integration 
formalism in Sec. lIIII and obtain the self-consistency (or- 
der parameter and number) equations. There we also dis- 
cuss the low energy scattering amplitude of a finite range 
interaction for all possible angular momentum channels, 
and relate the self-consistency equations to scattering pa- 
rameters. In Sec. lIVI we discuss the evolution from BCS 
to BEC superfluidity at zero temperature. There we ana- 
lyze the order parameter, chemical potential, quasiparti- 
cle excitation spectrum, momentum distribution, atomic 
compressibility and ground state energy as a function of 
scattering parameters. We also discuss gaussian fluctua- 
tions and low energy collective excitations at zero tem- 
perature in Sec.0 In Sec I VII we present the evolution of 
superfluidity from the BCS to the BEC regimes near the 
critical temperature. There we discuss the importance 
of gaussian fluctuations, and analyze the number of un- 
bound, scattering and bound fermions, critical temper- 
ature and chemical potential as a function of scattering 
parameters. In Sec. IVIII we derive TDGL equation and 
extract the GL coherence length and time. There, we re- 
cover the GL equation in the BCS and the CP equation 
in the BEC limit. A short summary of our conclusions is 
given in Sec. lVIIll Finally, we present in Appendices IX Al 
and IX BI the coefficients for the low frequency and long 
wavelength expansion of the action at zero and finite tem- 
peratures, respectively. 



II. GENERALIZED HAMILTONIAN 

The Hamiltonian for a dilute Fermi gas is given by 

k,Si k,k',q Si,S2,S3,S4 

6t^,,,(k,q)6,3,,Jk',q), (1) 
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where s„ labels the pseudo-spins corresponding to 
trapped hyperfine states and V is the volume. These 
states are represented by the creation operator aj^ , and 

^l,.2(k, q) = «lc+q/2,Si«-k+q/2,.2- ^^^^ = ^(^) " ^' 

where e(k) — fc^/(2A/) is the energy and /i is the chemical 
potential of fermions. 

The interaction term can be written in a separable 
form y/^3;,"/(k,k') = Tllf/^V{k,k'), where V/^f/^ is the 
spin and T^(k, k') is the spatial part, respectively. In 
the case of THS case, where s„ = (t,i), both pseudo- 
spin singlet and pseudo-spin triplet pairings are allowed. 
However, we concentrate on the pseudo-spin singlet THS 
state with rj3,s4 = Tll'^/J^i -8^3^,83^83 -S4- In addition, 
we discuss the SHS case (s„ =1), where only pseudo- 
spin triplet pairing is allowed, and the interaction is given 

by ^Tillt = rs?;sJ^si,s2^s3,s4^s4,T- this manuscript, we 
analyze THS singlet and SHS triplet cases for all allow- 
able angular momentum channels. THS triplet pairing 
is more involved due to the more complex nature of the 
vector order parameters, and therefore, we postpone this 
discussion for a future manuscript. 

The two fermion interaction can be expanded as 



V{k,k') = J d^rV{r)e 



i(k-k')T 



and should have the necessary symmetry under the Par- 
ity operation, where the transformation k —>■ — k or 
k' ^ -k' leads to V"(k,k') for singlet, and -V{k,k') 
for triplet pairing. Furthermore, l/(k, k') is invariant un- 
der the transformation (k, k') — > (— k, — k'), and V{k, k') 
reflects the Pauli exclusion principle. 

In order to obtain an approximate expression for the 
atomic interaction potential, we use the Fourier expan- 
sion of a plane wave in 3D 



An 



t'Mkr)Y,*{v)Y,^^,{k), 



(3) 



In the long wavelength limit (fc 0), one can show 
that the k dependence of this potential becomes exactly 
separable. In fact, for fcr <C 1 and k'r ^ 1, the asymp- 
totic expression of the spherical Bessel function for small 
arguments can be used, giving Ve{k,k') = Cek^k'^, with 
the coefficient Ci dependent on the particular choice of 
the real space potential. In the opposite limit, where 
kr ^ 1 and k'r ^ 1, the potential is not separable. In 
this case, Vi{k,k') mixes different k and k' , and shows 
an oscillatory behavior (which is dependent on the exact 
form of V{r)) with a decaying envelope that is propor- 
tional to l/{kk'). 

Under these circumstances, we choose to study a model 
potential that contains most of the features described 
above. One possibility is to retain only one of the £ terms 
in Eq. |0J), by assuming that the dominant contribution 
to the scattering process between fermionic atoms occurs 
in the ^th angular momentum channel. This assumption 
may be experimentally relevant since atom-atom dipole 
interactions split different angular momentum channels 
such that they may be tuned independently. Using the 
properties discussed above, we write 



Viik,k') = -XeTeik)Te{k'), 



(6) 



(2) where > is the interaction strength, and the function 



Tiik) 



{k/ko 



(1 



kVk^o) 



(7) 



describes the momentum dependence. Here, ko ^ Rq 
plays the role of the interaction range in real space and 
sets the scale at small and large momenta. In addition, 
the diluteness condition {uiRq <C 1) requires (ko/kp)^ » 
1 , where ng is the density of atoms and fcp is the Fermi 
momentum. This function reduces to T£{k) ~ k^ for 
small k, and behaves as Tg{k) ~ 1/k for large k, which 
guarantees the correct qualitative behavior expected for 
Ve{k, k') according to the analysis above. 



where ji{kr) is the spherical Bessel function of order £ 
and Yi^meO^) is the spherical harmonic of order (i,mi), 
in Eq.|2]to evaluate the matrix elements of the interaction 
potential in k-space 

V{k,k') = 47r ^ Veik, fc')>^£.n. (k)r,%^ (k')- (4) 

Here, J2e,m, = Efco Ei,=-£: and k denotes the an- 
gular dependence (0k,0k)- The {k,k') dependent coef- 
ficients Vi{k,k') are related to the real space potential 
V{r) through the relation 



III. FUNCTIONAL INTEGRAL FORMALISM 

In this section, we describe in detail the THS singlet 
case for even angular momentum states. A similar ap- 
proach for the SHS triplet case for odd angular momen- 
tum states can be found in Ref."^"', and therefore, we do 
not repeat the same analysis here. However, we point 
out the main differences between the two cases whenever 
it is necessary. 



A. THS Singlet Effective Action 



Ve{k,k')^47r drr^u{kr)3i{k'r)V{r). (5) 



The index £ labels angular momentum states in 3D, with 
^ = 0, 1, 2, ... corresponding to s,p, d, ... channels, respec- 
tively. 



In the imaginary-time functional integration formalism 
(h = k^ = 1, and (3 = ^/T), the partition function for 
the THS singlet case can be written as 



D{a\ 



a e 



'St 



(8) 
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with action 



Si 



dr 



Y.al^{T){dr)a^Ar) + Hi{T) 



(9) 



where the Hamiltonian for the ^th angular momentum 
channel is 

Hi{T) = ^eKk)at,,(r)ak,.(T) 



V 



Here, 6£,™,(q,T) = Ek rKfc)^^,™^ (k)ak+q/2,Tak-q/2a 
and ^^(k) — e(k) — ji^. We first introduce the Nambu 
spinor V'^(p) = (a], po.^p.i), where p = (k, iwj) denotes 
both momentum and fermionic Matsubara frequency 
Wj = (2j + l)7r//3, and use a Hubbard- Stratonovich trans- 
formation 



+fcT(g)<i>(q)-|-<ff(g)f,(g) 



(11) 



to decouple fermionic and bosonic degrees of freedom. 
Integration over the fermionic part [D{'tp\il>)] leads to 
the action 



J2 [Kd^)Sp,p' - Trln(G,//3)-i] , (12) 



p,p' 



where q = (q, ivj ) , with bosonic Matsubara frequency 
Vj = 2-11] / (3. Here, 



= $*(^)r,(^)a_ + $,(_q)r,(^)a+ 



-I- [iwjao - C«(k)CT3] 6p^ 



(13) 



is the inverse Nambu propagator, = 

Y.,ne^^,"ieiQ)^e,me{^^) is the bosonic field, and 

a± — ((Ti ± )/2 and at is the Pauli spin matrix. The 
bosonic field 



(14) 



has r-independent Ai^mt and r-dependent Ai^meiq) 
parts. 

Performing an expansion in 5'^'^ to quadratic order in 
A^,mf (g) leads to 

sr--sr + ^ E Ai„,(g)F-;„^„,(,)A,,„,,(,), 



(15) 



where the vector A] (g) is such that AJ (q) = 
[Al™,(g),A^,™,(-q)], and F-^^_,„,((j) are the matrix 

elements of the inverse fluctuation propagator matrix 
Fj^{q). Furthermore, is the saddle point action given 
by 



E[/?CKk)-Trln(Gf//3)-^ 



(16) 



and the saddle point inverse Nambu propagator is 

{Gfy^ =iwjao-Cii^)<J3 + A*,ik)a.+Ai{k)a+, (17) 
with saddle point order parameter 

A,(k) = r,(fc)^A,,™,r,,™,(k). (18) 

mi 

Notice that, A£(k) may involve several different m£ for a 
given angular momentum channel £. 

The matrix elements of the inverse fluctuation matrix 
Fj^{q) are given by 



CF"^ ,119 = 



^E(G?')"(f +^)(Gr)ii(|-fc) 



■(19) 



p 

r2(p)y,,„,(k)r;„,(k). ^^^^ 
^E(GJr)i2(f + ^)(G?^)i2(|-fc) 
r2(p)y,,„,,(k)r* ,(k). (20) 



Notice that while (Fjl^^^,Ji2{q) 
even under the transformations q - 



C^i.Li.m'^hiiq) are 
— q and ivj — > —ivj; 
(FZi,.,„pii(9) = (Ff"i^,,™P22(-(^) are even only under 
q — q, having no defined parity in ivj. 

The Gaussian action Eq. H15() leads to the thermody- 
namic potential $1?'^"'''' = + fi?"'^', where 



E 



^i,me I 



2 



AnV-^Xi 



^{C.(k)-i?,(k) 



ln[l + exp(-/3i?,(k))]}. 



= ^Ei^det[F7'(9)/(2/?)] 



(21) 
(22) 



are the saddle point and fiuctuation contributions, re- 
spectively. Here, 



£;,(k) = [ef(k) + |A,(k)p]- 



(23) 



is the quasiparticle energy spectrum. Having completed 
the presentation of the functional integral formalism, we 
discuss next the self-consistency equations for the order 
parameter and the chemical potential. 
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B. Self-consistency Equations 



The saddle point condition SSf^/SA*^ 



leads to 



the order parameter equation 



iTZ^ WIFTTX tanh , (24) 



V 



which can be expressed in terms of experimentally rele- 
vant parameters via the T-matrix approachSS.. 

The low energy two body scattering amplitude between 
a pair of fermions in the ^th angular momentum channel 
is given hy^ 



h{k) 



(25) 



where ri < and ag are the effective range and scattering 
parameter, respectively. Here rf has dimensions of L^^~^ 
and ag has dimensions of L^''+^, where L is the size of 
the system. The energy of the two body bound state is 
determined from the poles of fi{k ini), and is given 
by i?b,^ — — '«f/(2M). Bound states occur when oq > 
for £ = 0, and a^^^o^'^^o < for £ 7^ 0. Since ri < 0, 
bound states occur only when ai > for all £, in which 
case the binding energies are given by 



^'b.O — 



1 



Mai' 



1 



Matrt 



(26) 
(27) 



Notice that, only a single parameter (oq) is sufficient to 
describe the low energy two body problem for £ = 0, while 
two parameters (a^ , ri) are necessary to describe the same 
problem for £ 7^ 0. The point at which l/(fc^^+^a£) = 
corresponds to the threshold for the formation of a two 
body bound state in vacuum. Beyond this threshold, oq 
for £ = and ja^^o^^^^ol for £ 7^ are the size of the 
bound states. 

For any £, the two body scattering amplitude is related 
to the T-matrix via 



/,(fc) = -^r,[A:,fc;2e(k) + zO^ 
47r 



(28) 



where the T-matrix is given by 

. , , , 1 l/(k,k")T(k",k',S) 

T(k, k', E) = ViK k') + ^ E l,l2.(kV.O-. ■ 

Using the spherical harmonics expansion for both 
y(k, k') and T(k, k',i?) leads to two coupled equations, 

1 _ M 1 Tj{k) 



(29) 
(30) 



relating and k^ to and rp. Except for notational 
differences, notice that these relations are identical to 
previous results^. After performing momentum integra- 
tions we obtain 



kf^'+'a, 
1 



MkoXty^ 



MkoXi(j)i — AjTy/TT 



+ 2(£+l)0F 



(31) 
(32) 



where (t>i = r(£+l/2)/r(£+l) and (f)e = r(£-l/2)/r(£ + 
1). Here r(a;) is the Gamma function. Notice that. 



diverges and changes sign when MkgXi 



k^^^^L 

47rY^, which corresponds to the critical coupling for Fes- 
hbach resonances (the unitarity limit). 

300 



Mko?io 



150 








FIG. 1: Plot of original interaction strength MkoXo versus 
scattering parameter l/(fcoao). The inset shows MfcoAo versus 
l/(A;Fao) for ko « 200A:f. 



300 



Mko?ii 



150 






-0.25 







0.25 



0.5 



FIG. 2: Plot of original interaction strength MkoXi versus 
scattering parameter l/(fcoai). The inset shows MkgXi versus 
l/(fc|ai) for ko « 200A;f. 

In addition, the scattering parameter has a maximum 
value in the zero (Af — > 0) and a minimum value in the 
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infinite (A^ —f oo) coupling limits given respectively by 



J,2£+1 max _ 



t,2£+l min 



(ae > 0). 



{ae < 0), (33) 



(34) 



The first condition Eq. (|33|) (when 0) follows from 

Eq. (p?^ where rg^o < has to be satisfied for all possible 
ai^Q. However, there is no condition on ro for £ = 0, 
and koGQ^'^^ = in the BCS limit. The second condition 
Eq. H34(l (when — > oo) follows from Eq. 1)31(1 . which 
is valid for all possible £. The minimum ai for a finite 
range interaction is associated with the Pauli principle, 
which prevents two identical fermions to occupy the same 
state. Thus, while the scattering parameter can not be 
arbitrarily small for a finite range potential, it may go 
to zero as fco ^ oo. Furthermore, the binding energy is 
given by 



(35) 



when fco^+^af(^£ > 2(£ + 1)0F. 

In Fig. n we plot the original interaction strength 
MkoXo versus the scattering parameter kgao for the s- 
wave {£ = 0) channel. Notice that, fco|ao| in the 
BCS and fcoap ^ 1 in the BEC limit. A divergence in 
fcoflo corresponds to a s-wave Feshbach resonance occur- 
ring at MkoXo = An. 

In Fig. 13 we plot the original interaction strength 
MkoXi versus the scattering parameter /cqOi for the p- 
wave {£ — 1) channel. Notice that, fepjai | — 4 in the BCS 
and /cpai — » 2 in the BEC limit. A divergence in k^ai 
corresponds to a p-wave Feshbach resonance occurring at 
MfcoAi = Stt. 

Thus, the order parameter equation in terms of the 
scattering parameter is rewritten as 



MVA. 



167r2fc2^ 



E 



1 



t'Anh[/3Ei{k)/2] 



2e(k) 2Ee{k) 

,.r2(fc)r* (k)y,,„,(k). 



(36) 



This equation is valid for both THS pseudo-spin singlet 
and SHS pseudo-spin triplet states. However, there is 
one important difference between pseudo-spin singlet and 
pseudo-spin triplet states. For pseudo-spin singlet states, 
the order parameter is a scalar function of k, while it is 
a vector function for pseudo-spin triplet states discussed 
next. 

In general, the triplet order parameter can be written 
in the standard fori 



^40 



0,(k) 



-df{k) + idy{k) d|(k) 
df(k) 4(k) + i4(k) 



(37) 



where the vector dg{k) = [d^ (k) , d^ (k) , df (k)] is an 
odd function of k. Thefore, all up-up, down-down and 



up-down components may exist for a THS pseudo-spin 
triplet interaction. However, in the SHS pseudo-spin 
triplet case only the up-up or down-down component 
may exist leading to Ai(k) oc (Of)sisi(k). Thus, for 
the up-up case df{k) = and d^{k) = —id^{k), lead- 
ing to (i^(k) = df{k){l,i,0), which breaks time rever- 
sal symmetry, as expected from a fully spin polarized 
state. The corresponding down-down state has di{k) = 
(ij(k)(l, — i, 0). Furthermore, the simplified form of the 
SHS triplet order parameter allows a treatment similar 
to that of THS singlet states. However, it is important 
to mention that the THS triplet case can be investigated 
using our approach, but the treatment is more compli- 
cated. 

The order parameter equation has to be solved self- 
consistently with the number equation Ni = —dflt/dfii 
where fli is the full thermodynamic potential. In the 
approximations used. 



Nf,^N: 



gauss 



iVF + N, 



fluct 



(38) 



has two contributions. The saddle point contribution to 
the number equation is 



iVf = ^n,(k), 



(39) 



where ^^(k) is the momentum distribution given by 



neXk) 



Eiik) 



tanh 



f3Ee{k) 



(40) 



For the SHS triplet case, the summation over s is not 
present in N^^^ . The fluctuation contribution to the num- 
ber equation is 



detF.-i' 



(41) 



where ^ (g) is the inverse fluctuation matrix defined in 
Eq. CH) and (EOl). 

In the rest of the paper, we analyze analytically the 
superfluid properties at zero temperature (ground state) 
and near the critical temperatures for THS singlet (only 
even £) and SHS triplet (only odd ^) cases. In addition, 
we analyze numerically the s-wave (^ = 0) channel of 
THS singlet and p-wave (^ = 1) channel of SHS triplet 
cases, which are currently of intense theoretical and ex- 
perimental interest in ultracold Fermi atoms. 



IV. BCS TO BEC EVOLUTION AT T = 

At low temperatures, the saddle point self-consistent 
(order parameter and number) equations are sufficient 
to describe ground state properties in the weak coupling 
BCS and strong coupling BEC limitsi^. However, fluc- 
tuation corrections to the number equation may be im- 
portant in the intermediate regimeil. 
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Ground state properties {T = 0) are investigated by 
solving saddle point self-consistency (order parameter 
and number) equations to obtain A^.m^, and iii, which 
are discussed next. 



A. Order Parameter and Chemical Potential 

We discuss in this section A<i.m^ and fxg. In weak 
coupling, we first introduce a shell about the Fermi en- 
ergy < wd such that ep » wd ^ ^eC^p), in- 
side of which one may ignore the 3D density of states 
factor (-^Z e/ep) and outside of which one may ignore 
Af(k). While in sufficiently strong coupling, we use 
^£(k) 3> |A£(k)| to derive the analytic results discussed 
below. It is important to notice that, in strictly weak and 
strong coupling, the self-consistency equations Eq. H39() 
and l|36|) are decoupled, and play reversed roles. In weak 
(strong) coupling the order parameter equation deter- 
mines Ai^mt it^t) s-nd the number equation determines 

A^f (Ai^mt)- 

In weak coupling, the number equation Eq. H39|) leads 

to 



(42) 



for any i where ep = (2M) is the Fermi energy. In 
strong coupling, the order parameter equation Eq. 13t)|) 
leads to 



1 



(43) 
(44) 



where 
function 



= r(£- \/2)/T{i+l) and r(a;) is the Gamma 
This calculation requires that apfco > 1 for 
^ = and that fc^^+^a^0£ > (^ + 1)^^ for £ for 
the order parameter equation to have a solution with 
/if < in the strong coupling limit. In the BEG limit 
= ~kl/[2M{koao - if] for £ 0. Notice that, /iQ = 
-l/(2Mag) when k^ao > 1 [or |/io| < eo = kl/{2M)], 
and thus, we recover the contact potential (k^ oo) 
result. In the same spirit, to obtain the expressions in 
Eq. H43I) and H44I) . we assumed j/ifj ^ eo. Notice that, 
fii = i5b//2 in this limit for any £. 

On the other hand, the solution of the order parameter 
equation in the weak coupling limit is 



|Ao,o| = le^Trepexp 



fco 



exp 



TT /cp 

2-1 

2 fco 



2i-l 



2fcp|ao| 



,(45) 



2fcr+'l 



, (46) 



where h = 7r/4 and te>i = n2'^+^{2£ - 3)!!/^!. These 
expressions are valid only when the exponential terms 



are small. The solution of the number equation in the 
strong coupling limit is 



|A.,| ^8.p(g-) . 

E 2 64y^ _ 
|A£^0,mJ = -^^'^Fi'^F^o)^ 



(47) 
(48) 



to order iie/eo, where we assumed that ^^(k) 3> |A£(k)| 
for sufficiently strong couplings with j/i^ | <C eo- 




FIG. 3; Plot of reduced order parameter Ar — |Ao,o|/eF 
versus interaction strength l/(feFao) for fco « 200fcF. 




-2-1012 

l/lkpao) 

FIG. 4: Plot of reduced chemical potential fii — /to/ep (inset) 
versus interaction strength l/(fcFao) for fco ~ 200fcF. 

Next, we present numerical results for two particular 
states. First, we analyze the THS s-wave {£ = 0, me = 0) 
case, where Ao(k) = Ao,oro(fc)lo,o(k) with lo,o(k) — 
1/\/Att. Second, we discuss the SHS p-wave {£ — l^mg — 
0) case, where Ai(k) = Ai,ori(fc)Yi,o(k) with Yi,o(k) = 
a/3/ (47r) cos(0k)- In all numerical calculations, we choose 
fco ~ 200fcp to compare s-wave and p-wave cases. 
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In Figs, inland^ we show |Ao,o| and /io at T = 

for the s-wave case. Notice that the BCS to BEC evo- 
lution range in l/(fcFao) is of order 1. Furthermore, 
|Ao,o| grows continuously without saturation with in- 
creasing coupling, while /io changes from ep to £'b,o/2 
continuously and decreases as —l/{2Ma'^ for strong cou- 
plings. Thus, the evolution of |Ao,o| and /io as a func- 
tion of l/(fcFao) is smooth. For completeness, it is 
also possible to obtain analytical values of ao and Aq^o 
when the chemical potential vanishes. When /xq = 0, 
we obtain for |Ao,o| = 8eF[7r2\/^/r*(l/4)]i/3 « 3.73eF 
at l/(fcFao) = (27^3^eF/|Ao^o|)'/V[2r2(3/4)] « 0.554, 
which also agrees with the numerical results. Here T{x) 
is the Gamma function. 
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FIG. 5: Plots of reduced order parameter Ar = |Ai,o|/eF 
versus interaction strength l/(fc|ai) for fco ~ 200fcF. 
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FIG. 6: Plots of reduced chemical potential /ir = fii/ep 
(inset) versus interaction strength l/{kpai) for fco ~ 200fcF. 

In Figs. Inland IHl we show |Ai.o| and /ii at T = for 
the p-wave case. Notice that the BCS to BEC evolution 
range in l/(fc|ai) is of order k^/k-p. Furthermore, |Ai o| 
grows with increasing coupling but saturates for large 
l/(fcFai), while iii changes from ep to i?b,i/2 continu- 



ously and decreases as — l/(Mfcoai) for strong couplings. 
For completeness, we present the limiting expressions 



|Ai,o| = 24^eFexp 

Kf 



8 Trfeg TT 

'3 ifc^ ~ 2A;3|ai| 



5eF 



(49) 
(50) 



in the weak and strong coupling limits, respectively. 

The evolution of |Ai^o| and ni are qualitatively simi- 
lar to recent T = results for THS fermion^^ and SHS 
fermion-boson2& models. Due to the angular dependence 
of Ai(k), the quasiparticle excitation spectrum £^i(k) is 
gapless for /ii > 0, and fully gapped for /ii < 0. Further- 
more, both Ai.o and /ii are nonanalytic exactly when 
crosses the bottom of the fermion energy band /ii = 
at l/(A:|ai) « 0.48. The nonanalyticity does not occur 
in the first derivative of Ai o or /ii as it is the case in 
2D2^, but occurs in the second and higher derivatives. 
Thus, in the p-wave case, the BCS to BEC evolution is 
not a crossover, but a quantum phase transition occurs, 
as can be seen in the quasiparticle excitation spectrum 
to be discussed next. 



B. Quasiparticle Excitations 

The quasiparticle excitation spectrum Ei{}<l) = 
[^l(k) + |Af(k)p]^/^ is gapless at k-space regions where 
the conditions A^(k) = and e(k) = /i£ are both sat- 
isfied. Notice that the second condition is only satisfied 
in the BCS side fii > 0, and therefore, the excitation 
spectrum is always gapped in the BEC side {fie < 0). 

For £ = 0, the order parameter is isotropic in k-space 
without zeros (nodes) since it does not have any angu- 
lar dependence. Therefore, the quasiparticle excitation 
spectrum is fully gapped in both BCS (/iq > 0) and BEC 
(/io < 0) sides, since 



(51) 



min{£;o(k)} = |Ao(fc^J|, (/*o > 0), 
min{£;o(k)} = y'|Ao(0)|2+/i2, (^^ < Q). (52) 



Here, fc^^ = \/2MjIe. This implies that the evolution of 
the quasiparticle excitation spectrum from weak coupling 
BCS to strong coupling BEC regime is smooth when fiQ = 
for £ — pairing. 

In Fig. we show Eo{kx = 0,ky,kz) for an s-wave 
{£ = 0, = 0) superfluid when a) /io > (BCS side) 
for l/(fcpao) = —1, b) /io = (intermediate regime) 
for l/(fcpao) w 0.55, and c) /io < (BEC side) for 
l/(A:Fao) — 1. Notice that the quasiparticle excitation 
spectrum is gapped for all three cases. However, the sit- 
uation for £ ^ is very different as discussed next. 

For £ 0, the order parameter is anisotropic in k-space 
with zeros (nodes) since it has an angular dependence. 
Therefore, while the quasiparticle excitation spectrum is 
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FIG. 7: (Color online) Plots of quasiparticle excitation 
spectrum Eo{kx — 0,ky,kz) when a) jUo > (BCS side) 
for l/(A:Fao) — —1, h) fio — (intermediate regime) for 
l/(fcFao) « 0.55, and c) /io < (BEC side) for l/(feFao) = 1 
versus momentum ky/kp and k^/kp- 



gapless in the BCS {m^o > 0) side, it is fully gapped in 
the BEC ifie^Q < 0) side, since 

min{S,^o(k)} - 0, (/i^ > 0), (53) 
min{S,^o(k)} = (m£<0). (54) 



This implies that the evolution of quasiparticle excita- 
tion spectrum from weak coupling BCS to strong cou- 
pling BEC regime is not smooth for £ 7^ pairing hav- 
ing a nonanalytic behavior when fig-io — 0. This sig- 
nals a quantum phase transition from a gapless to a fully 
gapped state exactly when fj-e^o drops below the bottom 
of the energy band jie^o = 0. 

In Fig. IHl we show Ei{kx = 0,ky,kz) for a p-wave 
{£ = l,me = 0) superfluid when a) ^1 > (BCS side) 
for l/(fc|ai) = — 1, b) /ii = (intermediate regime) 
for l/(fc|ai) w 0.48, and c) /ii < (BEC side) for 
l/(A:pai) = 1. The quasiparticle excitation spectrum is 
gapless when Ai(k) oc kz/kp = and k"^ + ky + k1 = 
2M111 are both satisfied in certain regions of k-space. 
For kx = 0, these conditions are met only when kz = 
and ky = ±y/2Mfii for a given fii. Notice that, these 
points come closer as the interaction (/ii) increases (de- 
creases), and when — they become degenerate at 
k = 0. For /ii < 0, the second condition can not be sat- 
isfied, and thus, a gap opens in the excitation spectrum 
of quasiparticles as shown in Fig. |HJ:. 

The spectrum of quasiparticles plays an important role 
in the thermodynamic properties of the evolution from 
BCS to BEC regime at low temperatures. For £ ~ 0, 
thermodynamic quantities depend exponentially on T 
throughout the evolution. Thus, a smooth crossover oc- 
curs at /io = 0. However, for £ ^ 0, thermodynamic 
quantitites depend exponentially on T only in the BEC 
side, while they have a power law dependence on T in 
the BCS side. Thus, a non-analytic evolution occurs at 
fJ-e^o = 0. This can be seen best in the momentum dis- 
tribution which is discussed next. 



C. Momentum Distribution 

In this section, we analyze the momentum distribution 
neik) = [1 - ^i{k)/Ee(k)] /2 in the BCS (^^ > 0) and 
BEC sides {^e < 0), which reflect the gapless to gapped 
phase transition for nonzero angular momentum super- 
fiuids. 

In Fig. ini we show noikx — 0,ky,kz) for an s-wave 
{£ ~ 0, — 0) superfluid when a) /io > (BCS 
side) for l/(fcFao) = — li b) /io = (intermediate 
regime) for l/(fcFao) ~ 0.55, and c) /io < (BEC 
side) for l/(fcFao) = 1. As the interaction increases the 
Fermi sea with locus ^o(k) = is suppressed, and pairs 
of atoms with opposite momenta become more tightly 
bound. As a result, no(k) broadens in the BEC side 
since fermions with larger momentum participate in the 
formation of bound states. Notice that, the evolution is a 
crossover without any qualitative change. Furthermore, 
TiQ^kx, ky — 0, kz) and noikx, ky, kz — 0) can be trivially 
obtained from no{kx = 0,ky,kz), since no{kx,ky,kz) is 
symmetric in kx , ky and kz . 

In Fig. 1101 we show ni{kx — 0,ky,kz) for a p-wave 
{£ = l,me = 0) superfluid when a) /ii > (BCS side) 
for l/(fc|ai) — —1, b) /ii = (intermediate regime) 
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FIG. 8: (Color online) Plots of quasiparticle excitation 
spectrum E\[kx = 0,ky,kz) in a) /xi > (BCS side) 



for 1/ik-^ai) 



-1, b) = (intermediate regime) for 



l/(fci^ai) « 0.48, and c) /ii < (BEC side) for l/(k^ai) = 1 
versus momentum ky/kp and k^/kp- 




FIG. 9: (Color online) Contour plots of momentum dis- 
tribution no{kx = 0,ky,kz) when a) /xo > (BCS side) 
for l/(fcFao) — —1, h) fiQ = (intermediate regime) for 
l/(fcFao) « 0.55, and c) < (BEC side) for l/(A:Fao) = 1 
versus momentum ky/kp and kz/kp. 



for l/(fc|ai) « 0.48, and c) /ii < (BEC side) for 
l/(/c|ai) = 1. Notice that ni{kx = 0,ky,kz) is largest 
in the BCS side when k^/kp — 0, but it vanishes along 
kz/kp = for any ky/kp in the BEC side. As the in- 
teraction increases the Fermi sea with locus ^i(k) = 



is suppressed, and pairs of atoms with opposite mo- 
menta become more tightly bound. As a result, the 
large momentum distribution in the vicinity of k = 
splits into two peaks around finite k reflecting the p- 
wave symmetry of these tightly bound states. Further- 
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FIG. 10: (Color online) Contour plots of momentum distribu- 
tion ni(kx = 0, fcy, fcz) in a) > (BCS side) for l/(fc|ai) — 
— 1, b) ^1 = (intermediate regime) for l/(fc|ai) ~ 0.48, and 
c) Hi < (BEC side) for l/(fc|ai) — 1 versus momentum 
ky/kp and kz/kp. 



more, ni{kx,ky,kz = 0) = [1 - sgn($i(k))]/2 for any 
/ii, and ni{kx,ky — 0, fc^) is trivially obtained from 
ni{kx = Q,ky,kz), since ni(k) is symmetric in kx,ky. 
Here, sgn is the Sign function. 

Thus, ni(k) for the p-wave case has a major rearrange- 



ment in k-space with increasing interaction, in sharp con- 
trast to s-wave. This qualitative difference between p- 
wave and s-wave symmetries around k = explicitly 
shows a direct measurable consequence of the gapless to 
gapped quantum phase transition when fj,i = 0, since 
ni(k) depends explicitly on £'i(k). These quantum phase 
transitions are present in all nonzero angular momen- 
tum states, and can be further characterized through the 
atomic compressibility as discussed in the next section. 



D. Atomic Compressibility 

At finite temperatures, the isothermal atomic com- 
pressibility is defined by kJ{T) = -{dV / &P)t,nJV 
where V is the volume and V is the pressure of the gas. 
This can be rewritten as 



«F(r) 



1 



_ J_ f dNi 



where the partial derivative dNg/d^e at T 
by 



E 



|AKk)P 



, (55) 

T,V 

s is given 
(56) 



The expression above leads to Kg (0) — 2N{€f)/Nq in 
weak coupling BCS and kJ(0) = 2iV(eF)eF/(3|/io|A^o) in 
strong coupling BEC limit for £ = 0, where A^(eF) = 
MV/cF/(27r^) is the density of states per spin at the 
Fermi energy. Notice that Kq{0) decreases as ag in 
strong coupling since |^o| = l/(2A'^ao)- However, we 
only present the strong coupling results for higher an- 
gular momentum states since they exhibit an interesting 
dependence on ai and ko. In the case of THS pseudo- 
spin singlet, we obtain ^^^^(0) = AN{e-p)e-p<l)e./ {eo(f)tN'j) 
for I > 1, while in the case of SHS states we ob- 
«f(0) 



tam 



N{eF)eF/{^U^^I^Nj) for ^ = 1 and 



kJ^^(O) = 2N{tF)eF(t)il{f-G(t)iNf) for I > 1. Here - 
r(£-l/2)/r(£-hl) and ^£ = r(^ - 3/2)/r(£-h 1), where 
r(a;) is the Gamma function. Notice that decreases 
as y/oi for £ = 1 since — l(Affcoai) and k^~>i(0) is a 
constant for ^ > 1 in strong coupling. 

In Fig. ^2 we show the evolution of Ko'(O) for a s-wave 
{i = O^me = 0) superfluid from the BCS to the BEC 
regime. Ko{0) decreases continuously, and thus the evo- 
lution is a crossover (smooth) as can be seen in the inset 
where the numerical derivative of Kq (0) with respect to 
l/(A:Fao) is shown {dKg (0)/d[(A;Fao)^^]}- This decrease 
is associated with the increase of the gap of the excitation 
spectrum as a function of l/(A:Fao)- In this approxima- 
tion, the gas is incompressible [kq (0) — > 0] in the extreme 
BEC hmit. 

In Fig. ^1 we show the evolution of (0) for a p-wave 
{£ = l,m£ = 0) superfluid from the BCS to the BEC 
regime. Notice that, there is a change in qualitative be- 
havior when Hi =0 at l/(fcpai) 0.48 as can be seen in 
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FIG. 11: Plot of reduced isothermal atomic compressibil- 
ity Kr — k'o{0)/ko versus interaction strength l/(fcFao) for 
fco ~ 200fcp. The inset shows the numerical derivative of 
dKr/d[{kFao)~^] versus l/(fcFao). Here, kq is the weak cou- 
pling compressibility. 



the inset where the numerical derivative of k^(O) with re- 
spect to l/(fepai) is shown {dK^(0)/(i[(fcpai)~-^]}. Thus, 
the evolution from BCS to BEC is not a crossover, but a 
quantum phase transition occurs when /ii = 
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FIG. 12: Plot of reduced isothermal atomic compressibil- 
ity Kr ~ k.1{0)/ki versus interaction strength l/(fcpai) for 
fco ~ 200fcp. The inset shows the numerical derivative of 
dKj:/d[{kpai)~^] versus l/(fc|ai). Here ki is the weak cou- 
pling compressibility. 



limits. On the other hand, when fxg = 0, d N^'^ /d/ij is 
finite only for ^ = 0, and it diverges for ^ ^ 0. This di- 
vergence is logarithmic for i = 1, and of higher order for 
£ > 1. Thus, we conclude again that higher derivatives of 
iV^^ are nonanalytic when fJ-i^o = 0, and that a quantum 
phase transition occurs for £ 0. 

Theoretically, the calculation of the isothermal atomic 
compressibility kJ{T) is easier than the isentropic atomic 
compressibility nf{T). However, performing measure- 
ments of nf (T) may be simpler in cold Fermi gases, since 
the gas expansion upon release from the trap is expected 
to be nearly isentropic. Fortunately, nf{T) is related to 
kJ{T) via the thermodynamic relation 



Cf{T) 



(58) 



where kT{T) 



> K 

Cf{T) > Cf{T) 
{T « 0) the ratio Cf{T)/C^{T) 



^{T) since specific heat capacitites 
Furthermore, at low temperatures 



const, and therefore, 



K^(T « 0) oc (T « 0). Thus, we expect qualitatively 
similar behavior in both the isentropic and isothermal 
compressibilities at low temperatures (T « 0) . 

The measurement of the atomic compressibility could 
also be performed via an analysis of particle density fluc- 
tuations^Sii^. As it is well know from thermodynamics^, 
kJ{T) is connected to density fluctuations via the rela- 
tion 



(nf) - {my = TimyKj (T), 



(59) 



where {ng) is the average density of atoms. From the 
measurement of density fluctuations kJ{T) can be ex- 
tracted at any temperature T. 

It is important to emphasize that in this quantum 
phase transition at /i^^o = Oj the symmetry of the order 
parameter does not change as is typical in the Landau 
classification of phase transitions. However, a clear ther- 
modynamic signature occurs in derivatives of the com- 
pressibility suggesting that the phase transition is higher 
than second order according to Ehrenfest's classification. 
Therefore, if the symmetry of the order parameter does 
not change when fig changes sign, what is changing? To 
address this question, the topology of momentum space 
is discussed next. 



E. Topological quantum phase transitions 



The non-analytic behavior occurring when fie^o — 
can be understood from higher derivatives of Kg with re- 
spect to fj,i 



1 



T.V 



(57) 

For instance, the second derivative d'^N^^/dfi^ = 
SJ2k.s l^^(k)PC£(k)/[2i;|(k)] tends to zero in the weak 
[m « ep > 0) and strong (yU^ « -Eb,^/2 < 0) coupHng 



In what follows, we discuss the role of momentum space 
topolog}i22i^2ii^ in the non-analytic behavior of the ther- 
modynamic potential, when fii^o = 0. To investigate the 
role of topology, we make an immediate connection to the 
Lifshitz transition'*- in the context of ordinary metals at 
T = and high pressure. In the conventional Lifshitz 
transition, the Fermi surface e(k,V) — ep changes its 
topology as the pressure V is changed. For an isotropic 
pressure V, the deviation AV = V — Vc from the critical 
pressure Vc is proportional to A/i — jj, — Hc where /ic 
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is the critical chemical potential at the transition point. 
A typical example of the Lifshitz transition is the dis- 
ruption of a neck of the Fermi surface which leads to a 
non-analytic behavior of the number of states Af{fJ-) in- 
side the Fermi surface. In this case, A/'(/i) behaves as 
^(/ic)+S|/x — /LtcP^^ for fj, < fic, and as A{^c) for n > Hc, 
in the vicinity of Hc- Here, /C = {3/2)B\fj, — ^c\^/'^ /ri^ is 
the electronic compressibility, where ric — N^/V is the 
critical density of electrons at the transition point. No- 
tice that /C is nonanalytic, but it is not singular. The 
quantity that signals a phase transition in this case is 
not /C, but the thermopower Q, which is proportional to 
~d\n{n'^IC)/dfi, thus leading to Q oc — |A^|^^/^. In the 
conventional Lifshitz transition, the system lowers its en- 
ergy by A£ (X -|A/i|^/2 cx -\AVf^'^, and the transition 
is said to be of second and half order—. 

The topological transition proposed here is analogous 
to the Lifshitz transition in the sense that the sur- 
face in momentum space corresponding to Ei^^oi^) — 
Eg^oCk, Hi) — changes from a well defined set of k 
points for /i^^o > to a null set for fit^o < 0. Here, 
Ei^oCk, fig) plays the role of e(k, T') and fJ-i^o = = 
plays the role of the critical pressure Vc- 

For the Lifshitz transition in £ ^ superfluids, there 
is a non-analytic behavior in d'^Ni^^o/dfj.f, and thus in 
9k^q(0)/9/j,£. This behavior in k^q(O) is due to the 
collapse of all order parameter nodes at k = 0, which 
produce a gap in the excitation spectrum Ei^Q(k) and 
a massive discontinuous rearrangement of the momen- 
tum distribution n^^o(k) in the ground state as iJ,i^o 
f^l^o = 0. A direct topological analogy with the stan- 
dard Lifshitz transition can be made by noticing the col- 
lapse of locus of zero quasiparticle excitation energy at 
fJ-e^o = fJ-e^o in the excitation spectrum of the system. 
Generalized topological invariants can be invented along 
the lines of Ref—i^^ii^, however, we do not discuss them 
here. Instead, we analyze next the phase diagram at zero 
temperature. 



F. Phase diagram 

To have a full picture of the evolution from the BCS 
to the EEC limit at T = 0, it is important to analyze 
thermodynamic quantities at low temperatures. In par- 
ticular, it is important to determine the quantum crit- 
ical region (QCR) where a qualitative change occurs in 
quantities such as the specific heat, compressibility and 
spin susceptibility. Here, we do not discuss in detail the 
QCR, but we analyze the contributions from quasiparti- 
cle excitations to thermodynamic properties. However, 
the discussion can be extended to include collective exci- 
tationgSS (see Sec. 01. 

Next, we point out a major difference between i = 
and £ states in connection with the spectrum of the 
quasiparticle excitations (see Sec. II V 1?|) and their contri- 
bution to low temperature thermodynamics. 

For i = 0, quasiparticle excitations are gapped for 
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FIG. 13: The phase diagram of s-wave superfluids as a func- 
tion of l/(fcFao). 



all couplings, and therefore, thermodynamic quantitites 
such as atomic compressibility, specific heat and spin sus- 
ceptibility have an exponential dependence on the tem- 
perature and the minimum energy of quasiparticle excita- 
tions exp[-min{£;o(k)}/r]. Using Eqs. lEH and lO 
leads to ~ exp[-|Ao(fc^o)|/T] in the BCS side (^o > 0) 
and ~ exp[-A/|Ao(0) |2 t he BEC side (/iq < 0) 

as shown in Fig. El where fc^^ — \/2MjIe. Notice that, 
there is no qualitative change across /io = at small but 
finite temperatures. This indicates the absence of a QCR 
and confirms there is only a crossover for s-wave (£ = 0) 
superfluids at T = 0. 



(BCS side) 


^^1 


= (BEC side) 


n,>o 




Hi<0 


gapless 




gapped 


CiaT^i 




Ciaexp[-|Hi|/T] 




\ — 





(kin,) 



FIG. 14: The phase diagram of p-wave superfluids as a func- 
tion of l/(fc|ai). 



For £ ^ 0, quasiparticle excitations are gapless in the 
BCS side and are only gapped in the BEC side, and there- 
fore, while thermodynamic quantitites such as atomic 
compressibility, specific heat and spin susceptibility have 
power law dependences on the temperature ~ j^/3*^o [^i 
the BCS side, they have exponential dependences on the 
temperature and the minimum energy of quasiparticle 
excitations ~ exp[— min{£'^^o(k)}/r] in the BEC side. 
Here, (3i^o is a real number which depends on particular 
£ state. For £ = 1, using Eqs. (O and (|53)) leads to 
- T'^i in the BCS side (^ii > 0) and - exp(-|/^i|/r) 
in the BEC side {fii < 0) as shown in Fig. Notice 
the change in qualitative behavior across /ii = (as well 
as other £ states) at small but finite temperatures. 
This change occurs within the QCR and signals the ex- 
istence of a quantum phase transition (T = 0) for £ ^ 
superfluids. 
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G. Thermodynamic Potential 

Now, we discuss the thermodynamic potential fli at 
T = in the asmyptotic BCS and BEC limits. In the 
weak coupling limit, we obtain 



with 



5 



which is identical to the full thermodynamic potential 
^Ig. This indicates that fi^"'^' is negligible in the BCS 
limit. 

However, in the strong coupling limit, we obtain 



(61) 



Notice that, /iB.£ = — E'b,^ > is the Bosonic chemi- 
cal potential and N-^/ = Ni/2 is the number of bosons. 
To evaluate ^J.B,i^ it is necessary to find the first nonva- 
nishing correction for — E^^i. In the specific case of 
£ = 0, we obtain ^b,q = 4eFfcFao/(37r) — 47raB.o/A^B,o 
for the chemical potential and iT/ = —irNQaa/ [MV) = 
— -^b,oMb,o for the thermodynamic potential of the pairs. 
Here ob.o = 2ao and Mb,o = 2M is the scattering length 
and mass of the corresponding bosons. A better estimate 
for aB,o ~ 0.6ao can be found in the literatur o'^^i'^^i^°i^^ . 
The main reason for this difference is that our theory 
does not include possible intermediate (virtual) scatter- 
ing processes which renormalize a^fi- This is also the 
case when we analyze the collective modes in Sec. IV CI 
and the TDGL equation in Sec. IVII Bl 

Using ^g — {fiB,i + E\y.i)/2 and the thermodynamic 
relation fii — [dEi/ dNi)v , where Ei is the ground state 
energy, we obtain 



1 



(62) 



Notice that this expression is identical to the thermo- 
dynamic potential of bosons VIb,i — Eb/ — Nb/Hb/, 
where Eb,i is the ground state energy. Therefore, the 
fermionic thermodynamic potential in the strong cou- 
pling limit should lead to the thermodynamic potential 
of real bosons {Qe = ^b/)- Since = fi^^ -I- 51^"'^*, we 
conclude from this thermodynamic argument that 



1 



^i^^ — 2^B,fA^B,i 



(63) 



in the strong coupling limit. Therefore, Eg — NiE\^^g/2 = 
Eb,i, or Ei/Nt ~ m = {Ebz/Nb^i - Hb,i)I'^ which is 
consistent with quantum monte carlo calculations^. 

The fluctuation contribution to 17^"'^' comes from the 
zero point energy of the collective excitations, which is 
discussed next. 



V. GAUSSIAN FLUCTUATIONS 

The pole structure of Ff(q, TOj) determines the two- 
particle excitation spectrum of the superconducting state 



lO"*", and has to be taken into account 



to derive fi^"'^'. The matrix elements of F^(q, iwj) are 
F£,mf,m^(q7 ivj) for a given £. We focus here only on the 
zero temperature limit and analyse the collective phase 
modes. In this limit, we separate the diagonal matrix 
elements of F^^^ ^, {q) into even and odd contributions 



(60) with respect to ivj 



(■p-l \E 



E 



ii+i- + E+E^){E+ + E^] 



^ 2E+E4{iv,f-{E++E^Y] 
S V 

rKA:)y,,„.,(k)r;„,(k) - ^^^(64) 

-E- 



ii+i- + E+E^){tv,) 



2E+E-[{tVjY - {E++E-y] 

r2(fc)y,,„,,(k)r;„,(k). (65) 

The off-diagonal term is even in ivj, and it reduces to 

A+A^{E++E^) 



E 



2E+E_[{ivjy - {E++E^y 



r2(fc)y,,„,(k)r; .(k). 



(66) 



Here the labels ± denote that the corresponding variables 
are functions of k ± q/2. 

In order to obtain the collective mode spectrum, 
we express A£,m,(g) = T^,m,(<?)e'''*''"f = [pe^mAQ) + 
ide,me{q)]/V2 where n^mM)^ '&e.,mM)i Pi,mi{q) and 
(^e,mi{q) arc all real. Notice that the new fields 

P£.mA(l) = Ti,rne{q)cOs[dt,mtiQ)]^ and Og^raAl) = 

(q) sin[??^^m^(q)] can be regarded essentially as the 
amplitude and phase fields respectively, when "de^meil) 
is small. This change of basis can be described by the 
following unitary transformation 



1 



(:*0 









(67) 



From now on, we take as real without 

loss of generality. The diagonal elements of 
the fluctuation matrix in the rotated basis are 

i^Lnu,rn'^)ll + i^i.me.m'^)^^^ ^^^^ 
(F7^ ,) 



(F7^ ,)22 = (F7^ ,)fi ~ CF;^ ,)i2; and the off- 



diagonal elements are (F^ ,^^ ,^, )i2 = (F 

-1 NO 



')21 = 



i(F„ ^, with the q dependence being implicit. 



A. Collective (Goldstone) Modes 

The collective modes are determined by the poles of 
the propagator matrix F^ (g) for the pair fluctuation fields 
Ki^raiisi)-! which describe the Gaussian deviations about 
the saddle point order parameter. The poles of F^((7) 
are determined by the condition detF7^(q) = 0, which 
leads to 2(2^-1-1) collective (amplitude and phase) modes, 
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when the usual analytic continuation ivj —^'W + jO"*" is 
performed. Among them, there are 2i + 1 amplitude 
modes which we do not discuss here. 

The easiest way to get the phase collective modes is to 
integrate out the amplitude fields to obtain a phase-only 
effective action. Notice that, for £ ^ channels at any 
temperature, and for £ = channel at finite temperature, 
a well defined low frequency expansion is not possible for 
fig > due to Landau damping which causes the collec- 
tive modes to decay into the two quasiparticle continuum. 
A well defined expansion [collective mode dispersion w] 
must satisfy the following condition w <ti min{E^ + E^} . 
Thus, a zero temperature expansion is always possible 
when Landau damping is subdominant (underdamped 
regime). To obtain the long wavelength dispersions for 
the collective modes at T = 0, we expand the matrix 
elements of F7^ / to second order in Iql and w to get 



Dg^mt.m'W , 



Re 
iBp 



2 



qiqj 



(68) 

(69) 
(70) 



The expressions for the expansion coefficients are given 
in App. IX Al 

For £ = 0, the coefficients Cq^qq = Cofi_oSij and 

'3oo ~ Qofi.o^i.j are diagonal and isotropic in (i,j), 
and Po,o,o = vanishes. Here, Sij is the Kronecker delta. 
Thus, the collective mode is the isotropic Goldstone mode 
with dispersion 



Wo,o(q) 

Con 



Co,o|q|, 



^o,o,oQo,o,o 



^o,o,o-Ro,o,o + -Bq Q 



(71) 
(72) 



where Co,o is the speed of sound. Notice that the quasi- 
particle excitations are always fully gapped from weak 
to strong coupling, and thus, the Goldstone mode is not 
damped at T = for all couplings. 

For £ 0, the dispersion for collective modes is not 
easy to extract in general, and therefore, we consider the 
case when only one of the spherical harmonics Yg^mi (k) 
is dominant and characterizes the order parameter. In 
this case, Pe,mt,m,t = due to the order parameter equa- 
tion, and the collective mode is the anisotropic Goldstone 
mode with dispersion 



(73) 

1 
2 

(.74) 



Notice that the speed of sound has a tensor structure 
and is anisotropic. Furthermore, the quasiparticle exci- 
tations are gapless when fit=^o > 0, and thus, the Gold- 
stone mode is damped even at T = 0. However, Landau 
damping is subdominant and the real part of the pole 
dominates for small momenta. In addition, quasiparticle 
excitations are fully gapped when /i^/o < 0, and thus, 
the Goldstone mode is not damped. Therefore, the pole 
contribution to f^^^g* comes from the Goldstone mode for 
all couplings. In addition, there is also a branch cut rep- 
resenting the continuum of two particle scattering states, 
but the contribution from the Goldstone mode dominates 
at sufficiently low temperatures. 

It is also illustrative to analyze the eigenvectors of 
F^^(g) in the amplitude-phase representation corre- 
sponding to small We^rui (q) mode 







_ (^e,meiq) 





-W£,m,(q) 



1 



(75) 



Notice that, when Bg^mi,me the amplitude and phase 
modes are not mixed. 

Next, we discuss the dispersion of collective modes in 
the weak and strong coupling limits, where the expansion 
coefficients are analytically tractable for a fixed {£, mi) 
state. 



B. Weak coupling (BCS) regime 

The s-wave {£ = 0, mi = 0) weak coupling limit is char- 
acterized by the criteria fiQ > and //q ~ ep ^ l^o,o|- 
The expansion of the matrix elements to order |qp and 
w"^ is performed under the condition [w, |qp/(2A/)] ^ 
|Ao,o|- Analytic calculations are particularly simple in 
this case since all integrals for the coefficients needed 
to calculate the collective mode dispersions are peaked 
near the Fermi surface. We first introduce a shell about 
the Fermi energy |^o(k)| < iud such that ep ^ wb ^ 
Ao(kF), inside of which one may ignore the 3D density of 
states factor -y/e/ep and outside of which one may ignore 
Ao(k). In addition, we make use of the nearly perfect 
particle-hole symmetry, which forces integrals to vanish 
when their integrands are odd under the transformation 
fo(k) — ^o(k). For instance, the coefficient that couple 
phase and amplitude modes vanish (-Bo.0,0 = 0) in this 
limit. Thus, there is no mixing between phase and ampli- 
tude fields in weak coupling, as can be seen by inspection 
of the fiuctuation matrix Fo(<;). 

For £ = 0, the zeroth order coefficient is 

A^(ep) 



^0.0,0 — 



47r 



and the second order coefficients are 



■-"0,0,0 



Qofi,o _ N{eF)v 



3 

-Ro,o,o 



0,0,0 — 



36|Ao,o 
12IA, 



0,0| 



(76) 

(77) 
(78) 
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Here, vp = kp/M is the Fermi velocity and N{ep) ~ 
MVkp I (27r^) is the density of states per spin at the Fermi 
energy. 

In weak couphng, since S^,„^ <C Ae^rnt,mtRe 



the sound velocity is simplified 

[Q^'mf.mfM^,™*:™*] foi" ^- Using the coefficients 
above in Eq. O, for ^ = 0, we obtain 



to 



(79) 



which is the well known Anderson-Bogoliubov relation. 
For i ^ 0, the expansion coefficients require more de- 
tailed and lengthy analysis, and therefore, we do not dis- 
cuss here. On the other hand, the expansion coefficients 
can be calculated for any £ in the strong coupling BEG 
regime, which is discussed next. 



C. Strong coupling (BEC) regime 

The strong coupling limit is characterized by the crite- 
ria ^le < 0, l^ifl < eo = k^/{2M) and |Cf(k)| > |A^(k)|. 
The expansion of the matrix elements to order |qp and 
w'^ is performed under the condition [w, |qp/(2M)] <C 
The situation encountered here is very different 
from the weak coupling limit, because one can no longer 
invoke particle-hole symmetry to simplify the calculation 
of many of the coefficients appearing in the fluctuation 
matrix Fi(q). In particular, the coefficient B£ „it,m'^ 
indicates that the amplitude and phase fields are mixed. 



Furthermore, Pi 



, since this coefficient reduces 



to the order parameter equation in this limit. 
For £ = 0, the zeroth order coefficient is 



^0,0,0 



the first order coefficient is 



^|Ao,o|' 



(80) 



So,o,o = (81) 
and the second order coefficients are 

-Do, 0,0 = -Ro,o,o 



K 



(82) 
(83) 



where k — N{€F)/{32^/\jio\eF). 

Using the expressions above in Eq. ifT^ , we obtain the 
sound velocity 



the sound velocity is given by Cb,o = 47raB,o'^B,o/^^B o- 
Making the identification that the density of pairs is 
nB.Q = no/2, the mass of the pairs is A/b,o = 2M and 
that the Bose scattering length is ob.o = 2ao, it follows 
that Eq. is identical to the Bogohubov result Cb,o- 
Therefore, our result for the fermionic system represents 
in fact a weakly interacting Bose gas in the strong cou- 
pling limit. A better estimate for ob.o ~ 0.6ao can be 
found in the literature^SiiSiSi^. This is also the case 
when we construct the TDGL equation in Sec. IVII Bl 
For i 0, the zeroth order coefficient is 



A 



15(/)£K 

2eo0r 

the first order coefficient is 



|A^,™,p7f,{™,}, (85) 



(86) 



and the second order coefficients are 



l,m^,m' 



3k 

, ■Ome.m', 

SV'^oImiI 



(89) 



where k = N{e-p) / [2,2^/1^), = T{i - l/2)/r(^ + 1), 
= r(£-3/2)/r(£-hl) and 0^ = r(2^- 3/2)/r(2£ + 2). 
Here r(a;) is the Gamma function, and 71, is an angular 
averaged quantity defined in App. IX Bl 



In strong coupling, since 



> 



At^m(,meRi,va(,va(, the sound vclocity is simphfied 

to C'^;™. « [A^,^.,^M\U.raJBl^,,^,Y'^ for any L 
Using the expressions above in Eq. 174|l . for ^ 7^ 0, we 
obtain 



'^^ll,{mt}\^<i>l fcp 

4>i ko 



VF 



(90) 
(91) 



Therefore, the sound velocity is also very small and 
its smallness is controlled by the interaction range kg 
through the diluteness condition i.e. (/co/^f)^ ^ 1, for 
i 0. Notice that, the sound velocity is independent of 
the scattering parameter for £ 7^ 0. 

Now, we turn our attention to a numerical analysis 
of the phase collective modes during the evolution from 
weak coupling BCS to strong coupling BEC limits. 



Co,o — 



( |Ao,o| 
\32M\fio\TT 



VF 



37r 



(84) 



Notice that the sound velocity is very small and its small- 
ness is controlled by the scattering length oq. Further- 
more, in the theory of weakly interacting dilute Bose gas. 



D. Evolution from BCS to BEC regime 

We focus only on s-wave {£ — 0, — 0) and p-wave 
{£ = l,mi = 0) cases, since they may be the most rele- 
vant to current experiments involving ultracold atoms. 
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In Fig. El we show the evolution of Co,o as a func- 
tion of l/(fcFao) for s-wave case. The weak couphng 
Anderson-Bogohubov velocity Co,o — v-p/ V3 evolves con- 
tinuously to the strong coupling Bogoliubov velocity 
Co,o = i'F'\/^Fio/(37r). Notice that the sound velocity 
is a monotonically decreasing function of l/(fcFao)j and 
the evolution across /iq = is a crossover. 




FIG. 15: Plots of reduced Goldstone (sound) velocity 
(Co,o)r = Cofi/vF versus interaction strength l/(fcpao) for 
ko w 200A;f. 



In Fig. El we show the evolution of Cj'g as a function 
of l/(fcFai) for p-wave case. Notice that CJ'q is strongly 
anisotropic in weak coupling, since C^'q = C\'q « O.ddup 
and = \/3C[q « 0.79t;F, thus reflecting the order 
parameter symmetry. In addition, C^^q is isotropic in 

strong coupling, since Cj'g = UF^/Sfep/ (27rfco) ~ 0.049uf 
for fco ~ 200fcF, thus revealing the secondary role of the 
order parameter symmetry in this limit. The anisotropy 
is very small in the intermediate regime beyond /ii < 0. 
Notice also that, C^'q is a monotonically decreasing func- 
tion of l/(fcFai) in BCS side until jii = 0, where it sat- 
urates. However, CJ^'q = Cf'g is a nonmonotonic func- 
tion of l/(fcFai)^, and it also saturates beyond ^.i — 0. 
Therefore, the behavior of CI'q reflects the disapperance 
of nodes of the quasiparticle energy Ei (k) as changes 
sign. 

These collective excitations may contribute signifi- 
cantly to the thermodynamic potential, which is dis- 
cussed next. 



E. Corrections to fl^^ due to collective modes 

In this section, we analyze corrections to the saddle 
point thermodynamic potential fi^^ due to low energy 
collective excitations. The evaluation of bosonic Matsub- 
ara frequency sums in Eq. leads to O^"'^* —^ Qf^^, 
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FIG. 16: Plots of reduced Goldstone (sound) velocity 
(Cj'q )r — C^Q /vf (solid squares) and {Ci'Q)r ~ C^'q/vf (hol- 
low squares) versus interaction strength l/(fc|ai) for ko ~ 
200fcF. The inset zooms into the unitarity region. 



where 

f^r" = E V^q) + I In [1 - exp(-/3iy,(q))] } (92) 
q 

is the collective mode contribution to the thermodynamic 
potential. The prime on the summation indicates that 
a momentum cutoff is required since a long wavelength 
and low frequency approximation is used to derive the 
collective mode dispersion. Notice that the first term in 
Eq. (|92|) contributes to the ground state energy of the in- 
teracting Fermi system. This contribution is necessary to 
recover the ground state energy of the effective Bose sys- 
tem in the strong coupling limit as discussed in Sec II V Gl 
The corrections to the saddle point number equation 
Ni°^^ — —dn(°^^/dfi£ are due to the zero point mo- 
tion (A^l^) and thermal excitation {Nf'^) of the collective 
modes 



d 



q 



dWe{q) 
dfii 



riBiWiiq)]. 



(93) 



(94) 



Here nB{x) — l/[exp(/3a;) — 1] is the Bose distribution. 
For £ — 0, the last equation can be solved to obtain 
TV^ = -6(aCo,o/a/io)C(4)r4/(7r2C^ o), which vanishes at 
T — 0. Here C(a;) is the Zeta function. Similarly, A^^^g 
has a power law dependence on T, and therefore, van- 
ishes at T = since the collective modes are not excited. 
N^^ gives small contributions to the number equation in 
weak and strong couplings, but may lead to significant 
contributions in the intermediate regime for all £. The 
impact of N^^ on the order parameter and chemical po- 
tential in the intermediate regime may require a careful 
analysis of the full fluctuation contributionsH.. 
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Until now, we discussed the evolution of superfluidity 
from the BCS to the BEC regime at zero temperature. 
In the rest of the manuscript, we analyze the evolution 
of superfluidity from the BCS to the BEC limit at finite 
temperatures. 



VI. BCS TO BEC EVOLUTION NEAR T = T^,i 

In this section, we concentrate on physical properties 
near critical temperatures T = Tc^i. To calculate Tc^f, 
the self-consistency (order parameter and number) equa- 
tions have to be solved simultaneously. At T = Tc/, 
then A^^„ij = 0, and the saddle point order parameter 
equation Eq. reduces to 



MV 



1 tanh[e£(k)/(2rc,,)] 
2e(k) 2U^) 



(95) 

This expression is independent of mi since the interaction 
amplitude Xe depends only on £. Similarly, the saddle 
point number equation reduces to 



(96) 



where np{x) ~ l/[exp(/3a;) -I- 1] is the Fermi distribution. 
Notice that the summation over spins (s) is not present 
in the SHS case. It is important to emphasize that the 
inclusion of TV^"'^* around T( ~ Td is essential to produce 
the qualitatively correct physics with increasing coupling, 
as discussed next. 



A. Gaussian Fluctuations 

To evaluate the gaussian contribution to the thermo- 
dynamic potential, we sum over the ferniionic Matsubara 
frequencies in Eq. H22I) . and obtain the action 



^fluct _ P_ 



E AL„,('z)i,:™,,„^(9)A.,™^(<?), (97) 

.-1 Aq) = (F7I 
fluctuation propagator given by 



where L„ ,(q) — (F„ ^ ,)ii is the element of the 



1 - nF(^+) - ^^F(^-) 



r2(fc)y,,„,(k)r* ,(k). 



(98) 



This is the generalization of the t — case to ^ ^ 0, 
where £,± = ^^(k ± q/2). From S^"'^', we obtain the 
thermodynamic potential fi^^"^'' = + where 
Vt^ is the saddle point contribution with A£(k) = 0, and 



(99) 



is the fluctuation contribution. 

We evaluate the bosonic Matsubara frequency {ivj) 
sums by using contour integration, and determine the 
branch cut and pole terms. We use the phase shift 



(^^"'^*(q, w) = Arg[det ^^(q, iVj + iO+)] 
to replace detL^(g) in Eq. H99|) . leading to 



rfluct / 



E 

q 



dw 



-nB(«;)^"^*(q,H, 



(100) 



(101) 



^f"'=*(q,0) and 



where ip'l""\(\,w 
n-Q{x) — l/[exp(/3a;) — 1] is the Bose distribution. No- 
tice that, this equation is the generalization of the s-wave 
(^ = 0) caseiiii^ for £ 0. Furthermore, the phase shift 
can be written as (^"'^*(q, w) = ^^'^(q, w) -f ^)?''(q, u;), 
where 



(/?f(q,w) = c/?£(q, w)e(w - Wq), 



(102) 



is the branch cut (scattering) and ^^^(q, w) is the pole 
(bound state) contribution. Here, 8(x) is the Heaviside 
function, w* — Wq^ — 2/i£ with Wc^ = |qp/(4M) is the 
branch frequency and /i^ is the fermionic chemical po- 
tential. 

The branch cut (scattering) contribution to the ther- 
modynamic potential becomes 



q 



dw 



(103) 



For each q, the integrand is nonvanishing only for w > 
w* since (^|'^(q, w) = otherwise. Thus, the branch cut 
(scattering) contribution to the number equation Nf^ = 
—dflf^/d^e is given by 



m 



q 



dw 
n 



(104) 



where w = w + w^. 

When ai < 0, there are no bound states above Tc^i and 
N^'^ represents the correction due to scattering states. 
However, when > 0, there are bound states repre- 
sented by poles at w < w*. The pole (bound state) 
contribution to the number equation is 



bs 



= -E'^B[wHq)]?74q,W£(q)], 



(105) 



where We{ci) corresponds to the poles of ^{q) and 



71t[q,Wi{(l)] =Res{ 



ddetLJ^[q,We{cD]/d^ii 
detL7i[q,>V,(q)] 



} (106) 



is the residue. Heavy numerical calculations are nec- 
essary to flnd the poles as a function of q for all cou- 
plings. However, in sufficiently strong coupling, when 
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np{£^±) ^ 1 in Eq. (|98|l . the pole (bound state) contri- 
bution can be evaluated analytically by eliminating Xi in 
favor of the two body bound state energy £'b,£ in vacuum. 
Notice that, E^ij is related to the E'b/ obtained from the 
T-matrix approach, where multiple scattering events are 
included. However, they become identical in the dilute 
limit. _ 

A relation between and i?b,^ can be obtained by 
solving the Schroedinger equation for two fermions inter- 
acting via a pairing potential V{r). After Fourier trans- 
forming from real to momentum space, the Schroedinger 
equation for the pair wave function is 

2e(k)V'(k) + i ^ y(k, k')^(k') = ^bV(k). (107) 

k' 

Using the Fourier expansion of y(k, k') given in Eq. Q 
and choosing only the ^th angular momentum channel, 
we obtain 



V ^ 2e(k) - Ei,, 



(108) 



This expression relates -Bb/ < to in order to express 
Eq. (|lll|l in terms of binding energy Eh,e < 0. Notice 
that, this equation is similar to the order parameter equa- 
tion in the strong coupling limit {^i < and \^£\ ^ Tc^e), 
where 



1 



-E 



2e(k)-2/x/ 



(109) 



Therefore, fj,i — > Eb^i/2 as the coupling increases. 

Substitution of Eq. ((TU5|l in Eq. lfTUH|l yields the pole 
contribution which is given by yVe{(l) — Wq + E^^i — 2/i£, 
and the residue at this pole is ?7f[q, Wf(q)] — ~'2J2mi- 
Therefore, the bound state contribution to the phase shift 
in the sufficiently strong coupling limit is given by 

(^^"(q,iy) = 7re(w- Wq + MB,^), (HO) 
which leads to the bound state number equation 



(111) 



where fiB,e — 2//^ — i?b,f < is the chemical potential 
of the bosonic molecules. Notice that, Eq. (|lll|l is only 
valid for interaction strengths where fiB.e < 0. Thus, 
this expression can not be used over a region of coupling 
strengths where fiB,i is positive. 



B. Critical Temperature and Chemical Potential 



To obtain the evolution from BCS to BEC, the number 



and order parameter [Eq. (|^ ] equations have to be 
solved self-consistently for Tc/ and ii£. First, we analyze 
the number of unbound, scattering and bound fermions 
as a function of the scattering parameter for the s-wave 
{£ = 0) and p-wave {£ = 1) cases. 

In Fig. 1171 we plot different contributions to the num- 
ber equation as a function of l/{kpao) for the s-wave 
{£ = 0, = 0) case. Notice that, N^^ (iV^") dominates 
in weak (strong) coupling, while Nq'^ is the highest for 
intermediate couplings. Thus, all fermions are unbound 
in the strictly BCS limit (not shown in the figure), while 
all fermions are bound in the strictly BEC limit. 
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FIG. 17: Fractions of unbound Fg^ = Ng^/No, scattering 
^o'^ = N^VNo, bound ff' = A^o'VA'o fermions at T = Tc,o 
versus l/(fcFao) for ko ~ 200fcF. 

In Fig. ^1 we present plots of different contributions 
to the number equation as a function of l/(fc|ai) for 
the p-wave {£ — l,mi — 0) case. Notice also that, Nl^ 
(Ni^) dominates in weak (strong) coupling, while Nf"^ is 
the highest for intermediate couplings. Thus, again all 
fermions are unbound in the strictly BCS limit, while all 
fermions are bound in the strictly BEC limit. 

Therefore, the total fluctuation contribution Nf^ + Nf^ 
is negligible in weak coupling and N^^ is sufficient. How- 
ever, the inclusion of fluctuations is necessary for strong 
coupling to recover the physics of BEC. However, in the 
vicinity of the unitary limit [1 / (k^^'^^ ae) — > 0], our re- 
sults are not strictly applicable and should be regarded 
as qualitative. 

Next, we discuss the chemical potential and criti- 
cal temperature. In weak coupling, we introduce a 
shell about the Fermi energy |C^(k)| < wu, such that 
/i£ w ep ^ wd ^ Tc/. Then, in Eq. I|95|l . we set 
tanh[|^£(k)|/(2rc,£)] — 1 outside the shell and treat the 
integration within the shell as usual in the BCS theory. 
In strong coupling, we use that min[^f(k)] = Tc.e 
and set tanh[^f(k)/(2Tc.£)] = 1. Therefore, in strictly 
weak and strong coupling, the self-consistency equations 
are decoupled, and play reversed roles. In weak (strong) 
coupling the order parameter equation determines Tc^i 
{^i) and the number equation determines /i£ {Tc,i)- 
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FIG. 18: Fractions of unbound F^^^ = Nl'^/Ni, scattering 
psc ^ Nr/Ni, bound F}"' = N^'^Ni fermions at T = %,! 
versus l/(fc|ai) for ko « 200fcF. 
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In weak coupling, the number equation Ni leads 



to 



for any £. In strong coupling, the order parameter equa- 
tion leads to 



1 



Mo 



2Mal ' 



(114) 
(115) 



where (jye = T{i- l/2)/r(£+ 1) and r(a;) is the Gamma 
function. This calculation requires oofco > 1 for ^ = 0, 
and kQ^~^^ai(l)i > {£ + 1)^/Tr for £ 7^ for the order pa- 
rameter equation to have a solution with fi£ < 0. Fur- 
thermore, we assume \fig\ <C eo = k'^/{2M) to obtain 
Eqs. (|114|l and (|115|l . Notice that, = -Eb//2 in this 
limit. 

On the other hand, the solution of the order parameter 
equation in weak coupling is 



c,0 



— ep exp 

TT 



7- 



TT fcp 

2H 

2 fco 



TT 



ep exp 



2fcp|ao|_ 

TT 



(116) 
(117) 



where 7 « 0.577 is the Euler's constant, ti — tt/A and 
ti>i — 7r2^+^(2£ — 3)!!/i!!. These expressions are valid 
only when the exponential terms are small. In strong 
coupling, the number equation Ng, Nf^ leads to 



rpTHS _ 



2tt 



Mr 



E™. am 



0.218 



■ep, (118) 



where Mbj = 2M is the mass of the bosonic molecules. 
Here, nc = A:p/(37r^) is the density of fermions. For THS 



Fermi gases, we conclude that the BEG critical temper- 
ature of s-wave superfluids is the highest, and this tem- 
perature is reduced for higher angular momentum states. 
However, for SHS Fermi gases 



T; 



SHS 



2tt 



E™,C(3/2) 



0.137eF 



(119) 



where ng — fc|/(67r^) and C,{x) is the Zeta function. Here, 
the summation over mg is over degenerate spherical har- 
monics involved in the order parameter of the system, 
and can be at most V,^ — 2£ + 1. For SHS states, we 
conclude that the BEG critical temperature of p-wave su- 
perfluids is the highest, and this temperature is reduced 
for higher angular momentum states. 

For completeness, it is also possible to relate ae and 
Tc^i when chemical potential vanishes. When jig — 0, the 
solution of number equation Eq. 1112|l is highly nontrivial 
and it is difficult to find the value of the scattering pa- 
rameter a| at ii£ = 0. However, the critical temperature 
in terms of a| can be found analytically from Eq. (|95|l as 



(113) ( ^ 



1+^ 



{2-2-'+i)T{£+me+^y 



(120) 



to order Tc//eo, where eo = kQ/{2M) > Tc/. Notice 
that, this relation depends on kg only through a|. 

On the other hand, if temporal fluctuations are 
neglected, the solution for To,/ from the sad- 
dle point self-consistency equations is \Eh,i\ — 
2TQ,i In [3V^(To//ep)3/2/4] and fie = Eh,t/2 which leads 
to ' 



2 In 



(l^b,£|/eF 



(121) 



up to logarithmic accuracy. Therefore, Tqj grows with- 
out bound as the coupling increases. Within this cal- 
culation, the normal state for T > Tqj is described by 
unbound and nondegenerate fermions since A^(k) = 

and \fii>\/TQ^i - ln(|i?b.d/eF)^^^ > 1- Notice that the 
saddle point approximation becomes progressively worse 
with increasing coupling, since the formation of bound 
states is neglected. 

We emphasize that, there is no phase transition across 
Tq^i in strong coupling. However, this temperature is re- 
lated to the pair breaking or dissociation energy scale. 
To see this connection, we consider the chemical equilib- 
rium between nondegenerate unbound fermions (f) and 
bound pairs (b) such that 6 ^ / t +/ i for THS singlet 
states and 6 ^ / t +/ T for SHS triplet states. 

Notice that Tq^ is sufficiently high that the chem- 
ical potential of the bosons and the fermions satisfy 
I Mb I S> T and |/^t| 3> T at the temperature T of inter- 
est. Thus, both the unbound fermions (f) and molecules 
(b) can be treated as classical ideal gases. The equi- 
librium condition ^b = 2/if for these non-degenerate 
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gases may be written as T In {nb[27r/ (MbT)]^/^} — £^b,^ = 
2rin{nf[27r/(AffT)]3/2}, where rib K) is the boson 
(fermion) density, Mh {M[) is the boson (fermion) mass, 
and Eh.e is the binding energy of a bosonic molecule. 
The dissociation temperature above which some fraction 
of the bound pairs (molecules) are dissociated, is then 
found to be 



lii{\Ebj\/eF 



(122) 



where we dropped a few constants of order unity. There- 
fore, the logarithmic term is an entropic contribution 
which favors broken pairs and leads to a dissociation tem- 
perature considerably lower than the absolute value of 
binding energy |£'b/|- The analysis above gives insight 
into the logarithmic factor appearing in Eq. (|121|l since 
To^e 7dissoc,£/2. Thus, Tq./ is essentially the pair dis- 
sociation temperature of bound pairs (molecules), while 
Tc/ is the phase coherence temperature corresponding to 
BEC of bound pairs (bosonic molecules). 




FIG. 19: Plot of reduced critical temperature Tr — Tc.o/ep 
versus interaction strength l/(fcpao) at T — T^fi for ko ~ 
200A;f. 

In Fig.^1 we show T^fi for the s-wave {£ — 0, = 0) 
case. Notice that Tc^ grows from an exponential de- 
pendence in weak coupling to a constant in strong cou- 
pling with increasing interaction. Furthermore, the mean 
field Tq.o Siiid gaussian Tc^o are similar only in weak 
coupling, while Tq q increases without bound as Tq q ~ 
l/[(Mao)| ln(fcFao)|] m strong coupling. When /xq = 0, 
we also obtain analytically Tc^o/ep ~ 2.15/(A;Fao)^ from 
Eq. H120() . The hump in Tcfi around l/(fcFao) ~ 0.5 is 
similar to the those in RefJ^, and might be an artifact 
of the approximations used here. Thus, a more detailed 
self-consistent numerical analysis is needed to determine 
if this hump is real. 

In Fig. 1201 we show /iq for the s-wave case, where 
it changes from ep in weak coupling to £'b,o/2 = 
— l/(2Mao) in strong coupling. Notice that, at 
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FIG. 20: Plot reduced chemical potential /ir = fio/ep (inset) 
versus interaction strength l/(fcpao) at T = rc,o for ko ~ 
200fcF. 



is qualitatively similar to fj,o at T — 0, however, it is re- 
duced at Tc_o in weak coupling. Furthermore, /io changes 
sign at l/(fcFao) ~ 0.32. 
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FIG. 21: Plot of reduced critical temperature Tr = Tci/^f 
versus interaction strength l/(fcpai) at T = Tc,i for ko ~ 
200fcF. 

In Fig.|^ we show Td for the p-wave {£ — 1, mi — 0) 
case. Tc,i grows from an exponential dependence in weak 
coupling to a constant in strong coupling with increasing 
interaction. For completeness, we present the limiting 
expressions 



Tc,i = -epexp 

IT 



TcA - 



27r 



B.l 



ni 



C(3/2) 



= 0.137eF, 



(123) 



(124) 



in the weak and strong coupling limits, respectively. Fur- 
thermore, the mean field Tg i and gaussian i are sim- 
ilar only in weak coupling, while To.i increases with- 
out bound as Tq^i ~ l/[(Mfcoai)| ln(fcp/coai)|] in strong 
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coupling. When = 0, we also obtain analytically 
Tci/ep « 1.75/(fc^at)2/3 from Eq. (fHUI) . The hump in 
the intermediate regime is similar to the one found in 
fermion-boson model^. But to determine if this hump is 
real, it may be necessary to develop a fully self-consistent 
numerical calculation. 
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FIG. 22: Plot of reduced chemical potential /ir = fii/ep 
(inset) versus interaction strength l/(fc|.ai) at T = Tc,i for 
fco ~ 200fcF. 

In Fig. 1221 we show fii for the p-wave case, where 
it changes from ep in weak coupling to i?b,i/2 — 
— l/(Affcoai) in strong coupling. Notice that, at T^^i 
is both qualitatively and quantitatively similar to fj,i at 
T = 0. Furthermore, /xi changes sign at l/(A:|ai) « 0.02. 

For any given £, mean field and gaussian theories lead 
to similar results for Tc^i and Tq/ in the BCS regime, 
while they are very different in the EEC side. In the lat- 
ter case, Tq^i increases without bound, however, gaussian 
theory results in a constant critical temperature which 
coincides with the BEC temperature of bosons. Notice 
that the pseudogap region Tc^£ <T < Tg/ for £ — state 
is much larger than £ ^ states since Tqj^q grows faster 
than Tqj^q. Furthermore, similar humps in Tc^e around 
l/{k'^^^^ae) ~ are expected for any £ as shown for the 
s-wave and p-wave cases, however, whether these humps 
are physical or not may require a fully self-consistent nu- 
merical approach. 

As shown in this section, the frequency (temporal) de- 
pendence of fluctuations about the saddle point is cru- 
cial to describe adequately the bosonic degrees of free- 
dom that emerge with increasing coupling. In the next 
section, we derive the TDGL functional near Tcj to em- 
phasize further the importance of these fluctuations. 



VII. TDGL FUNCTIONAL NEAR Tc.i 

Our basic motivation here is to investigate the low fre- 
quency and long wavelength behavior of the order pa- 
rameter near Tc^i- To study the evolution of the time- 
dependent Ginzburg-Landau (TDGL) functional near 



Tc^i, we need to expand the effective action S^^ in 
Eq. l(T^ around Ai^^e = leading to 

Qcff _ c^P I Qga-uss 
D£ — ij£ -r Off 

{9.i,™*„} 

Alm,3 (93)A£,m,, {qi -92 + 93)- (125) 

Here, A£_,„^ (q) is the pairing fluctuation field. 

We first consider the static part of i^m^ m' il) ' ^^'^ 
expand it in powers of qi to get 



(q, 0) = Ui, 



E 



(126) 



Next, we consider the time-dependence of the 
TDGL equation, where it is necessary to expand 
L7^ ,(0,Wo) — /(0, 0) in powers of w after 

analytic continuation ivj — > w + iO+. 

In the X = (x, t) representation, the calculation above 
leads to the TDGL equation 



"mil ,"1^2 



E ^^Am,„ } (0) {x) 

d' 



id, 



A£,™,,(x)=0, (127) 



which is the generalization of the TDGL equation to 
higher momentum channels of THS singlet and SHS 
triplet states. Notice that, for THS triplet states, there 
may be extra gradient mixing textures and fourth or- 
der terms in the expansion^, which are not discussed 
here. All static and dynamic expansion coefficients are 
presented in Appendix IX Bl The condition det ai = 
with matrix elements ai^rnt-^,m(^ is the Thouless criterion, 
which leads to the order parameter equation. The coeffi- 
cient c^''', _ refiects a maior difference between £ = 

and £ ^ cases. While Cq g q — co,o,o<5j.j is isotropic 
in space, c^'^g is anisotropic, thus characterizing 

the anisotropy of the order parameter. The coefficient 
be,{mi^}(0) is positive and guarantees the stability of the 
theory. The coefficient di^rne^.me^ is a complex number. 
Its imaginary part reflects the decay of Cooper pairs into 
the two-particle continuum for ni > 0. However, for 
fi{ < 0, imaginary part of d^.mfj.mf^ vanishes and the 
behavior of the order parameter is propagating reflecting 
the presence of stable bound states. 

Next, we present the asymptotic forms of ai^rnij^,mi^', 



which are used 



to recover the usual Ginzburg-Landau (GL) equation 
for BCS superfluids in weak coupling and the Gross- 
Pitaevskii (CP) equation for a weakly interacting dilute 
Bose gas in strong coupling. 
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A. Weak coupling (BCS) regime 



B. Strong coupling (BEC) regime 



The weak coupling BCS regime is characterized by 
fj,i > and /i£ « ep 3> T^./. For any given £, we find 
the following values for the coefficients 



k7 In 



T 



^,7 

C 



7a 



1 TT 



4eF 8Tc. 



(128) 
(129) 
(130) 
(131) 



where = iV(eF)(eF/eo)7(47r) with iV(eF) = 

MVk-p / {"2,11^) is the density of states per spin at the 
Fermi energy. Here 5. 



is the Kronecker delta, and 
'^I'me me ^^'^ li,{me} E^rs angular averaged quantities 
defined in App. IX Bl Notice that the critical transition 
temperature is determined by det ae = 0. 

In the particular case, where only one of the spherical 
harmonics 1^,™^ (k) is dominant and characterizes the or- 
der parameter, we can rescale the pairing field as 



(x) 



'^«,{ma(o) 



(132) 



to obtain the conventional TDGL equation 



-Si - 



GL \ 



2yj2 I GL 



d_ 

dt 



(133) 



Here, e, = (Te,, - T)/T,j, with « 1, {U,me)i{T) 



GL length and r^,™^ = 



{&.rn,)ihi^ is the characteristic 



'ids 



^GL 



/ei is 



the characteristic GL time. 

In this limit, the GL coherence length is given by 
= [KU,™.C(3)/(4^^)]^/^(eF/rc,,), which 
makes {£,fmi)i ™uch larger than the interparticle spac- 
ing fcf^ There is a major difference between £ ^ Q and 



7^ pairings regarding (^ 



GL 

me / 



While Cq'-'o o = co,o,o^jj 
is isotropic, c^'^o.m,^,™,, = 4;m,^,m,/^i IS m general 
anisotropic in space (see App. IX Bp . Thus, (■^q^o)* 



isotropic and (C^o,mJj not. 

TpL^ ^ _i/(4gp) + 7r/(8rc,£) showin 



Furthermore, 

that the dynamics of '^J,y^i,{x) is overdamped reflecting 
the continuum of fcrmionic excitations into which a pair 
can decay. In addition, there is a small propagating term 
since there is no perfect particle-hole symmetry. As the 
coupling grows, the coefficient of the propagating term 
increases while that of the damping term vanishes for 
/i£ < 0. Thus, the mode is propagating in strong cou- 
pling reflecting the stability of the bound states against 
the two particle continuum. 



The strong coupling BEC regime is characterized by 
< and eo = kl/{2M) > |^f| > T^^g. For ^ = 0, we 
find the following coefficients 

ao,o.o = (2|a*o| - |-Bb,o|) , (134) 



&o,{o}(0) 



(^0,0,0 — 2kq, 
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where k% — iV(eF)/(64-y/eF|Aio|)- Similarly, for £ 0, we 
obtain 



24(1)1 (2\^e\ - \Eb,e\) S (138) 



&£#0,{mj„}(0) = 15je^{me„} 



~2^' 



'^e=i0,me^,me2 ~ ^^l't>l^mi.^,mi.^5i^j , 
dl^O.me-^ .me^ 



2K/i(j>iSmi^ ,me^ j 



where k|-^o 



(139) 

(140) 
(141) 

e = r{e- 

-1-2), where 

^E^O^me-^ ^me^ 

is isotropic in space for any £. Thus, the anisotropy of 
the order parameter plays a secondary role in the TDGL 
theory in this limit. 

In the particular case, where only one of the spherical 
harmonics Y^^mf (k) is dominant and characterizes the or- 
der parameter, we can rescale the pairing field as 



N{eF)/{6AyfH^). Here (j: 
l/2)/r(f + 1) and = r(2£ - 3/2)/r(2^ 
r(x) is the Gamma function. Notice that. 



^e,mei^) = \/'^i,rne,me-Ae.me{x), 



(142) 



to obtain the conventional Gross-Pitaevskii (CP) equa- 
tion 



d 



£,m); 



nrne = 



for a dilute gas 

/de 
tial, Mb/ 



2Mb I dt 

(143) 

of bosons. Here, fiB.i — 
2fii — Ei,j is the chemical poten- 
= 2M is the mass, 



Mdt^rai,mi/cY^ 



and Uii^rne = bi,{mi}{^)/d}„^^ „^^ is the repulsive inter- 
actions of the bosons. We obtain, f7o,o = 47rao/Af and 

C/<!#o,m, = 2^0'iT^^/^4>ai,{me}l{M4>lko) for ^ = and 
£ 7^ 0, respectively. Notice that the mass of the com- 
posite bosons is independent of the anisotropy and sym- 
metry of the order parameter for any given £. However, 
this is not the case for the repulsive interactions between 
bosons, which explicitly depends on £. 

For £ = 0, [/o,o = 47raB,o/AfB,o is directly proportional 
to the fermion (boson) scattering length ap (flB.o); where 
flB.o — 2ao is the boson-boson scattering lenth. A better 
estimate for a^fl ~ 0.6ao can be found in the litera- 
^^^.^48,49,50,51^ ^j^j^g foj. ^ _^ 0, Ue^me IS a Constant (inde- 
pendent of the scattering parameter a^) depending only 
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on the interaction range fcp and the particular {i,me) 
state. For a finite range potential, n-QjUi^m^ is small 
compared to ep, where n^,i — nf/2 is the density of 
bosons. In the £ = case nB,oC^,o/eF = 4fcFao/(37r) 
is much smaller than unity. For £ 7^ and even, 

nnjUt^raJe-F = 8O\/7i^0£7£,{mj}/(/)^(fcF/fco)- In the case 
of SHS states where (. ^ Q and odd, n-Q^iUi^mJ — 

^0\/TT(t>eli,{me} /4''e{f^F / ko)- The results for higher angu- 
lar momentum channels reflect the diluteness condition 
{kp/ko)^ « 1. 

To calculate {^fl^Ji in the strong coupling hmit, 
we need to know dfj,i/dT evaluated at Tc^i (see be- 
low). The temperature dependence of fig in the vicin- 
ity of Tc/ can be obtained by noticing that fiB,e — 
n{T)Ue,m,, where n{T) = nB,e [l - {T/T,^if/^] . This 
leads to kpi^fX)- = [^V(2MfcFf/f,mJ]'/' in the BEC 
regime. Using the asymptotic values of Ui,„l(^ we obtain 

kA& 

1/2 



[^/(Sfcpao)]'/' for £ = and kY{S.f^Q^^^)^ 



(/)2/(48OV^7,,{™j0,)J [ko/k^f/^ for I ^ 0. There- 
fore, (C^mj)i is also much larger than the interparticle 
spacing fcp^ in this limit, since fcpao for t — Q and 
fco ^ fcp for any 



C. Ginzburg-Landau Coherence length versus 
average Cooper pair size 

In the particular case, where only one of the spherical 
harmonics Y^^m^ (k) is dominant and characterizes the or- 
der parameter, we can define the GL coherence length as 



i.iLm,)i{T) 



/i2Mae, 



An expansion of 



the parameters ae^me,me and c^^^ in the vicinity of 
Tr f leads to 



large compared to interparticle spacing in both BCS and 
BEC limits, it should have a minimum near fig « 0. 
The prefactor of the GL coherence length must 

be compared with the average Cooper pair size S,^^" de- 
fined by 



pair \ 2 



(Z,(k)|r^|Z,(k)) 
(Z,(k)|Z,(k)) 
(Z,(k)|V^|Z,(k)) 
(Z,(k)|Z,(k)) 



{U7) 



where Zi{k.) — A£(k)/[2£^f (k)] is the zero temperature 
pair wave function. In the BCS limit, ^^^"^ is much larger 
than the interparticle distance fcp ^ since the Cooper pairs 
are weakly bound. Furthermore, for ^g < 0, we expect 
that is a decreasing function of interaction for any 
£, since Cooper pairs become more tightly bound as the 
interaction increases. Next, we compare (CFm^)* ^e^" 
for s-wave {£ — 0) and p-wave {£ = 1) states. 
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where the prefactor is the GL coherence length and given 
by 



/tGL ^2 



2MTc.g 



da. 



g,mi,m£ 



dT 



(145) 



The slope of the coefficient 
given by 



'.,mt,me with respect to T is 



da. 



'g,m£,rni 



E 



yg{ii) , ^^JLg ( 3^,(k) Xg{\^) 



2T2 dT \2T(g{k) Q{k) 



Here XgCk) and 3^f(k) are defined in App. IX Bl Notice 
that, while dfig/dT vanishes at Tc^g in weak coupling, it 
plays an important role in strong coupling. Furthermore, 
while (CFmi)j representing the phase coherence length is 



FIG. 23: Plots of GL coherence length kp^o.o (solid 
squares), and zero temperature Cooper pair size fcp^Q g'' (hol- 
low squares) versus interaction strength l/(fcFao) at T = Tc,o 
for fco ^ 200fcF. 

In Fig. 123 a comparison between (Co'o)i ^nd ^q^" 
is shown for s-wave {£ = 0,mg — 0). Co^" changes 
from kp^r" = [e'^/\/27r](eF/Tc,o) in the BCS limit 
to fcp^r'' = [eF/(2|Mo|)]^/^ = kFao/V2 in the BEC 
limit as the interaction increases. Here 7 ~ 0.577 
is the Euler's constant. Furthermore, when /jq = 0, 
we obtain kp^^"" = V7[T\l/4)/^^^yA w 1.29, 
where T{x) is the Gamma function. Notice that, ^q^" 
is continuous at /Lto = 0, and monotonically decreas- 
ing function of l/{k-pao) with a limiting value con- 
(146) trolled by ao in strong coupling. However, (^o^o)i i^ 
a non-monotonic function of l/(fcFao) having a mini- 
mum around l/{k-pao) « 0.32 (^0 = 0). It changes 
from A;f(Co°o)« - [7C(3)/(127r2)]i/2(,p/j; the BCS 

to fcF(Co:o)i = [7r/(8A:Fao)]^/^ in the BEC limit as the 
coupling increases, where (^(x) is the Zeta function. No- 
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tice that, (^^o)i grows as I/Vfcpoo in strong coupling 
limit. 
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FIG. 24; Plots of GL coherence length fcp^^o (solid 
squares), and zero temperature Cooper pair size fcp^i g"^ (hol- 
low squares) versus interaction strength l/(fcpai) atr = Tc,i 
for fco « 200fcF. 

In Fig. 131 a comparison between (^po)z ^^-'^ ^i"^"^ 
shown for p-wave {£ = l,m^ = 0). Notice that, (^f'^" is 
nonanalytic at fii = 0, and is a monotonically decreas- 
ing function of l/(fcpai) with a limiting value controlled 
by kp/ko in strong coupling. This nonanalytic behav- 
ior is associated with the change in £'i(k) from gapless 
(with hue nodes) in the BCS to fully gapped in the BEC 
side. However, (CPo)^ ^ non-monotonic function of 
l/(fc|ai) having a minimum around l/(fc|ai) w 0.02 
ifii = 0). It changes from kpi^^}^)^ = kF{^^}i)y = 
kF{^?fi)z/Vi = [7C(3)/(207r2)]i/2(ep/T^_^) the BCS 
to fcF(Cpo)» = [7rfco/(36fcF)]^/^ in the BEC limit as the 
couphng increases. Notice that, saturates in strong 
coupling limit reflecting the finite range of interactions. 

It is important to emphasize that (^p^f)z shown in 
Figs. H23|l and (|24|l is only qualitative in the intermediate 
regime around unitarity l/(fcp^^^a^) — since our theory 
is not strictly applicable in that region. 

VIII. CONCLUSIONS 

In this manuscript, we extended the s-wave {£ = 0) 
functional integral formalism to finite angular momen- 
tum £ including two hyperfine states (THS) pseudo- 
spin singlet and single hyperfine states (SHS) pseudo- 
spin triplet channels. We analyzed analytically super- 
fluid properties of a dilute Fermi gas in the ground state 
(T = 0) and near critical temperatures (T « Td) from 
weak coupling (BCS) to strong couphng (BEC) as a func- 
tion of scattering parameter (a^) for arbitrary i. How- 
ever, we presented numerical results only for THS s-wave 
and SHS p-wave symmetries which may be relevant for 



current experiments involving atomic Fermi gases. The 
main results of our paper are as follows. 

First, we analyzed the low energy scattering amplitude 
within a T-matrix approach. We found that bound states 
occur only when > for any £. The energy of the 
bound states i?b.f involves only the scattering parameter 
ao for ^ = 0. However, another parameter related to the 
interaction range 1/fco is necessary to characterize Eh,i 
for £ 0. Therefore, all superfluid properties for £ 7^ 
depend strongly on fcg and a^, while for £ — they depend 
strongly only on uq but weakly on fcg. 

Second, we discussed the order parameter, chemical 
potential, quasiparticle excitations, momentum distribu- 
tion, atomic compressibility, ground state energy, col- 
lective modes and average Cooper pair size at T = 0. 
There we showed that the evolution from BCS to BEC 
is just a crossover for £ = 0, while the same evolution 
for £ ^ exhibits a quantum phase transition character- 
ized by a gapless superfluid on the BCS side to a fully 
gapped superfluid on the BEC side. This transition is 
a many body effect and takes place exactly when chemi- 
cal potential fie^o crosses the bottom of the fermion band 
(/if = 0), and is best reflected as non- analytic behavior 
in the ground state atomic compressibility, momentum 
distribution and average Cooper pair size. 

Third, we discussed the critical temperature, chemical 
potential, and the number of unbound, scattering and 
bound fermions at T = Tc^e- We found that the criti- 
cal BEC temperature is the highest for £ = 0. We also 
derived the time-dependent Ginzburg-Landau functional 
(TDGL) near Tc/ and extracted the Ginzburg-Landau 
(GL) coherence length and time. We recovered the usual 
TDGL equation for BCS superfluids in the weak coupling 
limit, whereas in the strong coupling limit we recovered 
the Gross-Pitaevskii (CP) equation for a weakly inter- 
acting dilute Bose gas. The TDGL equation exhibits 
anisotropic coherence lengths for £ ^ which become 
isotropic only in the BEC limit, in sharp contrast to the 
£ = case, where the coherence length is isotropic for 
all couplings. Furthermore, the GL time is a complex 
number with a larger imaginary component for ni > 
reflecting the decay of Cooper pairs into the two particle 
continuum. However, for /i^ < the imaginary compo- 
nent vanishes and Cooper pairs become stable above T^j. 

In summary, the BCS to BEC evolution in higher an- 
gular momentum (£ 7^ 0) states exhibit quantum phase 
transitions and is much richer than in conventional £ ~ 
s-wave systems, where there is only a crossover. These 
£ 7^ states might be found not only in atomic Fermi 
gases, but also in nuclear (pairing in nuclei), astrophysics 
(neutron stars) and condensed matter (high- Tc and or- 
ganic superconductors) systems. 
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X. APPENDIX 
A. Expansion Coefflcients at T = 

From the rotated fluctuation matrix FJ^(q) expressed 
in the amphtude-phase basis, we can obtain the ex- 
pansion cocfhcicnts necessary to calculate the collective 
modes at T = 0. In the long wavelength (|q| 0), and 
low frequency limit (w 0) the condition 



{«',||}«min{2£;,(k)}, 



(148) 



is used. While there is no Landau damping and a well 
defined expansion is possible for i = case for all cou- 
plings, extra care is necessary for ^ 7^ when iig > Q 
since Landau damping is present. 

In all the expressions below we use the following sim- 
plifying notation Q = ^^^(k -t- q/2) /dqi, Ce'^ = ^^^(k -|- 
q/2)/(%9<7,), Al = 9A,(k + q/2)/% and A^'^' = 
d'^Ae{k + ci/2)/{dqidqj), which are evaluated at q = 0. 

The coefficients necessary to obtain the matrix element 
(F7"^ are 



A 



'-'mi 

- E J3r?m™.(k)i^,:„.^(k), (149) 



corresponding to the (q = 0, w = 0) term, 

- E Jj^{^>^'^^^k)[ef(k)-2A|(k)] 
+ 3A^'^"i?|(k)eKk)A,(k) + 5^]e^A2(k)e|(k) 

+ {ilA^+ijAl)Ae{-k) [2A2(k) - ^eA^)] } 

rf(fc)y,,„,(k)r;,„,(k), (i50) 

corresponding to the qiqj term, and 
eKk) 



E 



r,^(ft)r,,„,(k)y/,(k), (151) 



corresponding to the w'^ term. Here (5^^ is the Kro- 
necker delta. 

The coefficients necessary to obtain the matrix element 
(F7^ /)22 are 



P, 



'-'nil -^i^ 



+ A^'^"i?2(k)A,(k) 

+ 3^1^^AKk)+3A^A^,e£'(k) 

- mAl+ijADUmii^)} 

r2(fc)r,,™,(k)r;„,(k), (153) 

corresponding to the qiqj term, and 

Ri,m,,m', = E ^;^r2(fc)r,,™,(k)r,;„, (k), (i54) 

corresponding to the w'^ term. 

The coefficients necessary to obtain the matrix element 



E 



r2(fc)r,,„,(k)r* ,(k), (155) 



corresponding to the w term. 



B. Expansion Coefficients at T = Tc^e 

In this section, we perform a small q and ivj w+iO^ 
expansion near Tc/j where we assumed that the fluctua- 
tion field Ae_jni (x, t) is a slowly varying function of x and 
t. 

The zeroth order coefficient / (0, 0) is diagonal 



in m£ and m^, and is given by 



Aw 



--E 



A-^(k) 
2eHk) 



(156) 



where Ai(k) = tanh [/3^c(k)/2]. The second order coeffi- 
cient M&^LJ^^ ^1 {q,0)/{dqidqj) evaluated at q = is 
given by 



-E{ 



8^1 (k) 16eKk) 
/32pA',(k):y,(k) 



}r2(fc), (157) 



where 3^f(k) = sech^[/3^f(k)/2] and the angular average 

= / dkfc.fc,r,,™,(k)y,:^,(k). (158) 

Here, dk = sin(0k)rf^k'^'/>k, — sin(0k) cos(^k), ky = 



V--^r,^(fc)y,,„,(k)r,:„,(k){152) and = cos{9u). in general, a\ 



2Ei{k) 

corresponding to the (q = 0, w = 0) term. 



is a fourth order tensor for fixed £. However, in the 
particular case where only one of the spherical harmon- 
ics Yi^meO^) is dominant and characterizes the order pa- 



rameter, a„ 



is diagonal in mi 



and m'l. In this case, we use Gaunt coefficients^ to 
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show that -J^^ is also diagonal in i and j leading 



to a J , 



•^fjTnjjTnf "mf.mj^ij'- 



The coefficient of fourth order term is approximated at 
q„ = 0, and given by 



j(o) 



An 



E 



r^(fc),(i59) 



4^1 (k) seKk). 

where the angular average 

lUrni„} = J C^kYf,m^i W^^tm,^ (k) 

Y,^,n,,{k)Y,*^^Jk). (160) 

we expand 



To extract the time-dependence 



,(0,0) in 



powers of w after the analytic continuation ivj 
We use the relation {x ± iO+)~^ = P(l/x) q= i7r(5(a;), 
where V is the principal value and 6{x) is the Delta 
function, to obtain 



Qi 



4CKk) 



(161) 



Keeping only the first order terms in w leads to 

Ql,mf:,m'^{w + iO+) = -dl^rne,Tn'^W + ■■; whcrC 



47r 



E 

L k 



+ ^!^iV(,p) &2(M.)e(A*,) 
8 V 



(162) 



is also diagonal in mi and m^. Here iV(eF) = 
MVfcF/(27r^) is the density of states per spin at the Fermi 
energy, r^(a;) = {e^x^) / {e^ + xY^^ is the interaction sym- 
metry in terms of energy and Q{x) is the Heaviside func- 
tion. 
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